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O ■ We demonstrate that weakly coupled, large A^, d-dimensional SU (N) gauge theories on a 
class of compact spatial manifolds (including S'^~^x time) undergo deconfinement phase 
JlLi' transitions at temperatures proportional to the inverse length scale of the manifold in 
^ . question. The low temperature phase has a free energy of order one, and is characterized 
by a stringy (Hagedorn) growth in its density of states. The high temperature phase 
has a free energy of order A^^. These phases are separated either by a single first order 
transition that generically occurs below the Hagedorn temperature or by two continuous 
phase transitions, the first of which occurs at the Hagedorn temperature. These phase 
transitions could perhaps be continuously connected to the usual flat space deconfinement 
transition in the case of confining gauge theories, and to the Hawking-Page nucleation 
of AdS^ black holes in the case of the A/" = 4 supersymmetric Yang-Mills theory. We 
suggest that deconfinement transitions may generally be interpreted in terms of black hole 
formation in a dual string theory. Our analysis proceeds by first reducing the Yang-Mills 
partition function to a (0 -|- 0)-dimensional integral over a unitary matrix U, which is the 
holonomy (Wilson loop) of the gauge field around the thermal time circle in Euclidean 
space; deconfinement transitions are large A^ transitions in this matrix integral. 
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1. Introduction 

The thermodynamics of large SU{N) gauge theories is interesting for at least two 
different reasons. On the one hand, these theories are believed to share many qualitative 
features with the finite non-Abelian gauge theories which are relevant to real-world QCD 
(or GUTs), including confinement at zero temperature and a deconfinement transition as 
the temperature is increased. On the other hand, large gauge theories are believed 
to be dual to weakly coupled string theories, which display Hagedorn behavior (reviewed 
below) associated with singularities in various thermodynamic quantities. Thus, by 
studying large gauge theory thermodynamics, one may hope to achieve a better un- 
derstanding of both deconfinement transitions in gauge theory and Hagedorn behavior in 
string theory and, possibly, some relationship between the two. 

Unfortunately, asymptotically free gauge theories in Minkowski space are strongly 
coupled except at very high temperatures, so direct analysis of the thermodynamic behavior 
usually relies on numerical or lattice techniques. However, by placing the theory on a 
compact space, one obtains a tunable dimensionless parameter RAqcd (where R is the 
size of the compact space and Aqcd is the dynamically generated scale of the gauge 
theory0) which gives back the Minkowski space theory when it is large, but gives a weakly 
coupled theory when it is smallll. In this weak coupling limit the thermodynamics can 
be studied at all temperatures in perturbation theory, and this analysis is the goal of the 
present paper. Specifically, we investigate the thermodynamics of dimensional U{N) (or 
SU{N)) gauge theories with arbitrary adjoint- valued matter fields! on compact spatial 
manifolds, such as S'^~^, in the limit of weak coupling. This theory is simple enough to 
analyze exactly in the large A^ limit; we will find that it has a surprisingly rich structure 
with many of the expected features of strongly coupled gauge theory thermodynamics, 
including stringy Hagedorn behavior and a deconfinement transition. 

We should start by emphasizing that even at very weak coupling, SU{N) Yang-Mills 
theory on a compact manifold behaves very differently from (A^^ — 1) copies of the U{1) 

^ We use here, and throughout the paper, notations which are appropriate for four dimensional 
gauge theories. For lower dimensional gauge theories, Aqcd should be replaced by a scale which 
is an appropriate power of the 't Hooft coupling constant Oym^^ ■ 

^ A similar dimensionless parameter appears also in conformal gauge theories like the d = 4 
Af = 4 supersymmetric Yang-Mills theory, which can also be studied at weak coupling. 

^ Our methods can be easily generalized to include also matter fields in other representations. 
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theory. The excitations of the U{1) theory are arbitrary numbers of photons; in contrast 
(for example) a single SU{N) gluon is not an allowed excitation on a compact manifold, 
as gluons source lines of color electric flux. Due to Gauss' law, the only allowed excitations 
of an SU{N) theory on any compact manifold are combinations of gluons (and any other 
charged particles in the theory) that are grouped into SU{N) singlets. 

As we will see below, the projection onto the singlet sector introduces effective inter- 
actions between the gluons. It will turn out that these interactions are negligible at high 
energies, but they dominate the dynamics at low energies. As a consequence, in the N ^ oo 
limit, the gauge theories under study in this paper have at least two distinct phases. The 
lowest temperature phase (dominated by states with E/N"^ <^1/R) has a strongly stringy 
flavor; it is characterized (see for the free theory) by a Hagedorn-like density of states 
p{E) that grows exponentially with energy, p{E) oc e^/'^« with Th = C/R, where C IS a 
constant depending on the matter content of the theory and on the shape of the compact- 
ification manifold. On the other hand, the high temperature phase (dominated by states 
with E/N"^ ^ 1/R) behaves qualitatively like a gas of free particles; the free energy F 
takes the form F = N^f{T), where /(T) oc for T > In some cases these two 

phases may be separated by a more mysterious intermediate temperature phase. 

The free gauge theory undergoes precisely one phase transition as a function of the 
temperature; this transition is of first order and it occurs precisely at the Hagedorn tem- 
perature Th- As far as we know, this result was first found by ||^ in a specific case, and 
was derived more generally in a beautiful paper by Sundborg For finite A^, this phase 
transition is smoothed out, as it must be for any theory with a finite number of degrees of 
freedom at finite volume. 

Turning on a small 't Hooft coupling A = Qym-^ qualitatively alters the behavior of 
the free gauge theory. It turns out that the interacting theory displays one of two possible 
behaviors, depending on the sign of a coefficient determined by two-loop and three-loop 
vacuum diagrams. When varying the temperature, this theory either undergoes a single 
first order transition below the Hagedorn temperature, or it undergoes two continuous phase 
transitions, one at the Hagedorn temperature and the second one slightly above it. Note 
that, in the first case, the phase transition shields the Hagedorn spectrum of the theory in 

^ Sundborg's results do not seem to have received wide attention; in particular, we learned 
of his paper only after we had independently rederived all of its results, which have substantial 
overlap with section 3 and parts of section 5 of our paper. 
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the sense that there is no temperature at which string excitations of arbitrarily high energy 
dominate the partition function. In the second case the first of the two phase transitions is 
truly Hagedorn-like; the singularities in the neighborhood of the phase transition directly 
encode the high-energy density of states of stringy oscillators.ll 

We derive the results described in the previous paragraphs by reducing the thermal 
partition function to an integral over a single unitary U{N) (or SU{N)) matrix 



where U = e*^" (a is the zero mode oi Aq on Ai x , A4 is the spatial manifold in 
question, and (3 = l/T). The mode a is the lightest modei of the gauge theory on 
X S^] Seff(U), the quantum effective action for this light mode, may be computed by 
integrating out all other fields in the theory. This procedure (integrating out the heavy 
modes) may be explicitly performed in perturbation theory, generating an expression for 
Seff{U) as a power series in the 't Hooft coupling A. The lowest term in this power 
series, corresponding to the free gauge theory, is particularly easy to compute, either by 
evaluating one-loop vacuum diagrams or by explicitly counting gauge-invariant states in 
free Yang-Mills theory using a projection onto the singlet sector of the theory.0 For the 
free theory we obtain an expression (of the form (|1 . 1|) ) for Z{(3) which is exact, even at 
finite A^, and applies to any gauge theory with any matter content. The computation of 
higher order terms is straightforward in principle but tedious in practice. 

Having computed Seff we proceed to evaluate ( |1.1| ) in the usual manner, by changing 
the integration variable to the eigenvalues e*^ of U. For the adjoint theories we consider, 
the resulting effective action in the free limit is simply the sum over a pairwise potential 
between the eigenvalues with a temperature-independent repulsive term and an attractive 
term that increases from zero to infinite strength as the temperature increases from zero 
to infinity. In the large limit, the integral is determined by a saddle point characterized 

^ This is consistent with the general arguments of that a large N second order decon- 

finement transition must always be accompanied by a Hagedorn-like spectrum and occur at the 
Hagedorn temperature. 

^ We assume that there are no additional zero modes for any of the fields. This is not always 
true, since, for instance, the gauge field may have additional zero modes when the compactification 
manifold is not simply connected. We will not discuss such examples here. 

^ This second method was already used to derive similar expressions in [^J^. 
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by a density of eigenvalues p{6) that minimizes the effective action. At sufficiently low 
temperatures, the repulsive term dominates and p{6) is constant on the circle. At high 
enough temperatures, the attractive term forces the eigenvalues to bunch together and 
p{9) vanishes outside a narrow interval on the circle 0. In the free theory, these two 
regimes are separated by a first order transition at the Hagedorn temperature in which the 
eigenvalue distribution jumps discontinuously from the uniform distribution to a sinusoidal 
distribution. However, it turns out that depending on the details of the interactions in 
Seff, the nature of this transition may be modified at arbitrarily weak non-zero coupling. 

Based on the general form of 5'e//, we show that two classes of behavior are possible 
at weak coupling, as noted earlier. In the first case, we have a single first order phase 
transition at which the uniform low temperature distribution jumps discontinuously to a 
non-uniform distribution. This is similar to the behavior in the free theory, but the transi- 
tion happens strictly below the Hagedorn temperature for non-zero coupling. The second 
possibility is that the eigenvalue distribution evolves continuously as a function of temper- 
ature, changing from a uniform to a non-uniform (but nowhere vanishing) distribution at 
a first critical temperature, and then developing a gap (on which the distribution vanishes) 
at a second critical temperature. Consequently, the corresponding theory undergoes two 
phase transitions as a function of temperature; the first of these is a second order transition 
at the Hagedorn temperature, while the second is a third order transition closely related 
to the Gross- Witten phase transition of two dimensional lattice gauge theories. 

Which of these two possibilities is realized depends on the sign of a particular coeffi- 
cient in Sef f which appears to depend on the details of the field content and interactions 
of the theory in question. We are currently involved in a computation of this coefficient 
for the d = 4 pure Yang-Mills theory and for the d = 4 JV = 4 supersymmetric Yang-Mills 
(SYM) theory; we hope to report the result of this computation soon [ p!0| . 

The phase transition at (or near) the Hagedorn temperature has a natural interpreta- 
tion in the gauge theoryi; it is simply a deconfining transition. The low-temperature phase 
can be thought of as a gas of singlet 'glueballs', with the free energy scaling as A^^, while 
the high temperature phase, with the free energy scaling as A^, is a plasma of gluons (and 
other particles). Indeed, the traditional deconfinement order parameter, the Polyakov loop 
in the fundamental representation ^(tr(Pe^^^")) ~ ^(tr(t/)) , is a good order parameter 

^ As we have described above, in some situations the theory undergoes two phase transitions 
upon raising the temperature. The discussion in this paragraph appUes to the first of these. 
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for the phase transitions described in this paper. It is zero in the confining phase, and 
non-zero in the deconfining phase!. 

The phase transitions we discuss in this paper could potentiaUy be continuously related 
to several other interesting phase transitions - at least, the values of all the order param- 
eters are consistent with such a continuous relation. In the case of the 3 -|- 1-dimensional 
N = A SU {N) SYM theory compactified on a sphere, our analysis applies to the deconfine- 
ment transition at weak 't Hooft coupling; the deconfinement transition at strong coupling 
is related by the AdS/CFT correspondence [|12| to the Hawking-Page transition [|13| of 
gravitational theories on asymptotically anti-de Sitter (AdS) spaces. For asymptotically 
free gauge theories, our weakly coupled results apply when the scale R of the compactifica- 
tion manifold is much smaller than the strong coupling scale 1/Aqcd- However, when we 
take RAqcd ^ 1? our phase transitions could turn into the infinite volume deconfinement 
transition. 

We propose a set of phase diagrams for gauge theories as a function of coupling (A or 
RAqcd, respectively) and temperature that interpolate between our weakly coupled results 
and the known strong coupling behavior. We are led to speculate on a dual interpretation 
of deconfinement transitions; in particular, we conjecture that deconfinement transitions 
are always associated with black hole formation in a dual string theory. We argue that the 
mysterious intermediate temperature phase associated with second order deconfinement 
transitions (if such a phase is realized in some theory) would be dual to a string theory in 
a background dominated by a strange new type of stable black hole. 

The paper is organized as follows. We begin in §2 with a review of some of the relevant 
background concerning Hagedorn behavior in string theory, deconfinement transitions, 
and the relation between them. In §3, we show that large free gauge theories have a 
Hagedorn-like behavior of their spectrum and provide a first derivation of the matrix model 
expression for the thermal partition function of compactified gauge theories by explicitly 
summing over gauge-invariant states. In §4, we derive the same expression by using a 
path integral formalism. In §5, we describe the solution of this matrix model, exhibit the 
phase transition, and both compute and discuss the relevant order parameters. In §6, we 



^ Strictly speaking, the Polyakov loop is zero in both phases as flux conservation forbids placing 

a single quark on a compact space. In §5.7 we actually define our order parameter more carefully; 

it is either ^(tr(C7)) in an infinitesimal deformation of the gauge theory corresponding to an 

infinitesimal Higgsing of the theory, in the limit that the deformation goes to zero, or we can use 
^(tr([/)tr(^t)). 
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generalize the path integral derivation of ( p. . 1|) to interacting gauge theories and determine 
the general form of Seff in perturbation theory. Utilizing these results, we discuss the 
different possibilities for the phase structure of the theory at weak coupling. In §7, we 
suggest extrapolations of our results to strong coupling and propose a set of possible phase 
diagrams for gauge theories as a function of coupling and temperature. We also speculate 
on the interpretation of deconfinement transitions in terms of black hole formation in a 
dual string theory. We end in §8 with a summary of our results and a discussion of some 
possible future directions. Three appendices contain useful technical results. 



2. Background 

2.1. Hagedorn behavior in string theory 

A single free closed string has an infinite number of vibration modes, each of which 
may be excited to arbitrary level. Each vibrational state of the string corresponds to a 
distinct particle species in space-time. In fact (for strings in Minkowski space), the number 
of particle species grows exponentially with their mass. 

To be specific, consider a free type II superstring whose worldsheet is the direct sum of 
the free SCFT on 1R'^~^'^ and a compact unitary SCFT C with central charge c = 10 — d. 
Let N{M) denote the number of particle species in this theory with mass less than M. 
For much bigger than the string tension l/2na', 

dN Xe'^'^^^^ 

- — = J h subleadinq, (2.1) 

dM ^ ^ 

where X is a constant. In the zero coupling limit, each particle species is described by 
its own free quantum field theory. The finite temperature partition function Z of such a 
system of free quantum fields is easily computed; the contribution of highly excited string 
states (particle species of large mass, a'M'^ ^ 1) is given by0 



d-l 



M(27rv/2^-i) 

ln(Z(T)) = K'T"^ I dM ^ (2.2) 

M 2 



Note that d in ( p. I]) is the number of noncompact directions on the space IR''^^'^ x C' . 
It is possible that ( |2.1| ) (and the other formulas of this section) apply to more general string 
backgrounds (for instance with a warped space-time geometry) upon replacing d by deff, the 
effective number of noncompact space-time dimensions in such spaces. In particular, defj = 1 on 
a space like AdSs in global coordinates on which particle propagation is gapped. 
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for some constant K' . Note that Z{T) diverges for T larger than the Hagedorn temperature 
Th = „ AsTai 
develops a singularity 



= ——j=. As T approaches this critical temperature from below, the partition function 



ln(Z(r)) oc { (2.3) 

L (oi ) 2 d even 

where = - T. 

It is useful to understand the origin of the divergence in ( p.2|) from an alternate point of 
view. Recall that, for any system, the thermal partition function may be computed by the 
Euclidean partition function with the time direction compactified on a thermal circle S\ 
of circumference /3 = 1/T (bosons have periodic boundary conditions and fermions have 
anti-periodic boundary conditions around a thermal circle). For free string theory, this 
means that ^(T) in (|2.2| ) can be computed by the string theory torus partition function on 
the Euclidean space- timelR'^-^x^S^xC". If we now take the time direction to be in IR 
we can reinterpret ^(T) as the one- loop contribution to the vacuum energy of string theory 
on the 9 dimensional space IR'^"^ ^S\^C'. The one-loop contribution to the vacuum 
energy of any theory is determined completely by its spectrum; in particular it diverges 
if the spectrum includes a tachyon. Indeed, superstrings winding a thermal circle an odd 



number of times have a tachyonic mode when the size of the circle is small enough [|I4[ ; this 
is because modular invariance forces the imposition of the opposite GSO projection, one 
that projects in the identity operator, on states of odd winding. Specifically, the ground 
state of a superstring that winds once around the thermal circle has mass 

™^w = — \ — ^ — ^ - 1 I , (2.4) 

and becomes tachyonic for T > Th, providing an alternate explanation for the divergence 
of Z{T) for T > Th. 

For the theory with strictly vanishing string coupling Qs the partition function is simply 
ill-defined at T > Th- However, for any finite Qs, the perturbative analysis of the spectrum 



Note that the expectation value of the energy density at T = Th is convergent for d > 3 but 
diverges when d < 3. Note also that states of very high energy contribute significantly just below 
the Hagedorn temperature. Consequently, the properties of the theory at fixed energy E\foi' 
(microcanonical ensemble) may differ significantly from the properties of the theory at any fixed 
temperature (canonical ensemble). 
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breaks down at sufficiently liigli energies (of order l/gf), so tlie beliavior as T — > Th is 
difficult to analyze directly (recall that thermodynamics near the Hagedorn temperature 
receives non- negligible contributions from states with arbitrarily high energy). It has been 
speculated that at any finite coupling, no matter how small, the divergence of the free 



Z[j3) above the Hagedorn temperature may be replaced by a phase transition to a new high 
temperature phase. According to this suggestion, the string mode W winding the thermal 
circle is the order parameter for this transition [1T5|JT^ . The phase transition occurs when 



this mode condenses. This could occur either below the Hagedorn temperature (which is 
the temperature at which this mode becomes massless), in which case the transition is of 
first order, or at the Hagedorn temperature, in which case the transition is second order. 

Unfortunately, in general we know very little about this conjectured high temperature 
'phase' of string theory. In it was argued that, were such a phase transition to occur 



in flat space, it must be of first order. However, thermodynamics and phase transitions 
in gravitational systems in flat space are at best approximate notions (see [0 for a nice 
discussion of this issue). The high temperature 'phase' in flat space seems likely to be ill 
defined; in this phase energy densities are C'(p-) and would (since the Jeans instability is 
triggered at these energy densities) appear to involve black holes. As the density of states 
of Schwarzschild black holes grows even faster than exponential, the existence of such a 
phase seems problematic. 

The picture is much clearer for type IIB string theory on AdS^ x S^. At small 
curvatures, the spectrum of single string states in this theory exhibits a Hagedorn growth, 
with Hagedorn temperature Th = and d^ff = 1 (see §7.2). However, this theory 

(with temperature conjugate to global time) undergoes a first order phase transitionlll 
well below its Hagedorn temperature [1T3|,pT[]. As in flat space, energy densities in the 



high temperature phase are O(p-), and the high temperature phase is dominated by a 
big black hole sitting at the center of AdS space. In Euclidean space this black hole 
may indeed be thought of as a condensate of winding modeJH, so the Atick-Witten order 



As the space is effectively compact, tlie phase transition is sharp only in the limit Qs — > 0. 

Recall that the Euclidean time cycle is contractible in any black hole. Thus the 'time winding 
number' symmetry, present in global AdS, has been spontaneously broken in the black hole phase, 
implying that the formation of the black hole must involve the condensation of winding modes. 
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parameter analysis seems to apply, at least qualitatively, to this situation.0 Note that the 
entropy of big black holes grows relatively slowly with energy in AdS space (5* oc E^), and 
thermodynamics is well defined at all temperatures. 

In this paper we demonstrate that the thermodynamics of weakly coupled large A^, 
SU{N) gauge theories at finite volume (for instance gauge theories on a sphere) has some 
striking similarities to string thermodynamics on the weakly curved AdS space described 
above. The theories we study all undergo a phase transition as a function of temperature. 
The low temperature phase has a Hagedorn growth in the density of states with d^ff = 1. 
Finally, the high temperature phase in Euclidean space may be thought of as a condensate 
of winding modes. 

2.2. Deconfinement transitions in gauge theories 

In this paper we will study the thermodynamics of gauge theories on compact man- 
ifolds of size R. On taking R to infinity our theories reduce to gauge theories on fiat 



space, whose thermodynamics has been studied extensively (see [|T6[ for a review). In this 
subsection we review the thermal behavior of gauge theories on IR'^"^'^. We start with a 
discussion of confining gauge theories, and at the end of this section we briefiy mention 
the situation with the d = 4 JV = 4 supersymmetric Yang-Mills theory. 

It is widely believed that pure SU {N) gauge theories in d = 4 confine at zero tem- 
perature (see fl^ and for reviews); the low energy spectrum of the theory may be 
thought of as a Fock space of interacting glueballs.lll Consequently, it is expected that in 
such theories the low temperature (T ^ Aqcd) dynamics may be understood in terms of 
a sparsely populated thermal bath of glueballs. On the other hand, asymptotic freedom 
permits a reliable computation of the T ^ Aqcd behavior of such theorieslll, revealing 



Note, however, that it is not possible to use this analysis to predict (in this case) the order of 
the transition, as the Atick-Witten analysis, applied to string theories with a mass gap, could be 
consistent with either a first order or a second order phase transition, depending on the coupling 
constants of the theory (see appendix C). 

Similar expectations hold for gauge theories in lower dimensions, and for four dimensional 
gauge theories with a sufficiently small amount of charged matter fields; for QCD, the SU{3) 
gauge theory describing the strong interactions, they are experimentally verified. For simplicity 
of terminology we will restrict the discussion of this subsection to pure gauge theories in d = 4. 

In infinite volume this statement is far from trivial due to infrared divergences which afi^ect 
finite-temperature perturbative computations. These problems do not appear at finite volume so 
they will not be relevant for us. 
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that this high temperature phase may be understood as a weakly coupled gas of gluons 
1^,|19[. The high and low temperature phases appear qualitatively different; we will now 



review the argument showing that they are distinguished by an order parameter, and so 
are separated by a phase transition at some intermediate temperature. 

It is a defining property of a confining theory (with no fields in the fundamental 
representation of SU{N)) that a single external particle in the fundamental representation 
(a "quark") can be inserted into this theory only at an infinite cost in energy. Heuristically, 
such a quark forms one end of a QCD string which is infinitely long because it has nowhere 
else to end. In the low temperature confining phase this long string also has infinite 
positive free energy. Let Fq{T) represent the Yang-Mills free energy at temperature T in 
the presence of an external quark. It follows from the form of the coupling of an external 
quark to the gauge fields that e""^'*^^^/^ = (P), where V = ;^trPexp[— ^ A] is the so 
called Polyakov loop (sometimes called the Polyakov-Susskind loop |T^,^), the trace of a 
Wilson loop around the compactified Euclidean thermal time circle. Thus, {V) = in the 
low temperature confining phase. 

On the other hand, asymptotic freedom permits reliable computations that establish 
{V) 7^ at high enough temperatures. Thus, {V) constitutes an order parameter that 
sharply distinguishes the low temperature confining phase from the high temperature de- 
confined phase. From a low temperature point of view, the transition that separates these 
two phases is associated with the condensation of flux tubes whose effective free energy is 
driven negative at high enough temperatures (when the energy of these strings is overcome 
by the entropy of their vibrations). 

In the large limit (with fixed 't Hooft coupling, or equivalently fixed Aqcd) 
the deconfinement phase transition has yet another order parameter [^|. The confined 
phase is dominated by gauge-invariant bound states and so its free energy F{T) scales 
like A^^ at large A^. On the other hand, the deconfined phase is described by free glu- 
ons, and consequently its free energy scales as A^^ at large A^. Thus, in the large A^ 
limit, limjv^oo F{T)/N'^ constitutes a second order parameter for deconfinement; like the 
Polyakov loop this new order parameter vanishes in the confined phase but is finite in the 
deconfined phase. 

Next, we turn to a discussion of the order of the deconfinement phase transition. To 
study the phase transition one may generate an effective action for the Polyakov loop V 
by integrating all other fields out of the path integral that generates the Yang-Mills free 
energy. P is a complex scalar field on IR"^; note also that P — > is a symmetry of 
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the SU{N) Yang-Mills path integral (it is generated by gauge transformations that 
are single valued only up to an element of the center of the gauge group "En; such gauge 
transformations act trivially on the dynamical fields of the pure Yang- Mills theory), which 
is spontaneously broken in the deconfinement phase transition. 

The effective action for V completely determines the nature and properties of the 
deconfinement transition. If this transition turns out to be first order then its properties 
depend on the details of the theory under consideration. However, if the deconfinement 
transition is of second order then the phase transition point has universal behavior; it must 
be described by a ^Tv-invariant fixed point of a complex scalar field on IR"^ [p3| , p4| . Only 
a small number of fixed points with the required properties are known. For N = 2 the 
conditions above uniquely pick out the Wilson Fischer fixed point (the d = 3 Ising model). 
For = 3 no attractive fixed point with the symmetries listed above is known. For > 4 
the only known fixed point is U{1) invariant (all operators of the form (p'^^ for integer k 
and > 4 are irrelevant at this fixed point, hence the enhanced symmetry). 

The arguments of the paragraph above apply to all confining SU{N) gauge theories 
with adjoint matter. In the rest of this subsection, however, we will describe what is known 
about the deconfinement transition of the pure gauge theory on IR'^'^. Since this phase 
transition is a strong coupling phenomenon, we cannot study it in perturbation theory, and 
most knowledge about it comes from lattice simulations. Lattice simulations indicate that 
the SU{2) deconfinement phase transition is second order and is indeed in the Ising model 
universality class. The SU{3) deconfinement phase transition is first order, in agreement 
with the predictions above. However, the SU{N) deconfinement phase transitions for 
> 4 (at least for A^ = 4,5,6) also appear to be first order; in particular the U{1) 
invariant [d = 3 XY fixed point) appears not to be attained in pure Yang -Mills0 ( see, for 
instance, [|2^, but see also |l26|]). 

As the final topic in this subsection we review the thermodynamics of the A/" = 4 
SYM theory in fiat space. We focus on the limit N ^ oo with A = Qym^ held fixed. As 
this theory is conformal (and from large A^ counting) the free energy density F must take 



the form F = —^N f{X)T . Perturbative computations establish that, at small A, 



However, it is entirely possible that this beautiful fixed point describes the phase transition of 
some other confining gauge theory - say, for instance, TV = 1 supersymmetric Yang-Mihs theory. 
The critical exponents of this fixed point imply that the singular piece of the free Energy scales like 
F oc {5T)^'^^ (where ST = T — Tc) upon approaching the phase transition, implying deff ~ 5.16, 
see ( |2.3| ). It would be fascinating to find a stringy interpretation or confirmation of this result. 
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/(A) = 1 — while computations using the AdS/CFT correspondence show that, at 
large A, /(A) = | + It is generally assumed (see, e.g., [|9l) that /(A) is analytic 

on the positive real line (i.e. that one encounters no phase transitions in extrapolating 
the A/" = 4 SYM coupling from zero to infinity; more about this in §7). In summary, the 
thermodynamics of the A/" = 4 SYM theory in fiat space is qualitatively rather boring. 
The same theory has much more interesting thermodynamics on a sphere as was pointed 
out in [011; we will discuss this issue in great detail below. 



2.3. Hagedorn versus deconfinement 

't Hooft has argued 0] that it should be possible to recast any gauge theory, in the 
limit of large with fixed Qym-^^ ^ string theory with string coupling Qs oc If we 

accept this conjecture, then it is reasonable to expect that the gauge theory deconfinement 
transition has a dual description as a stringy thermal transition of the sort described in 
52.10. Indeed, stringy and gauge thermodynamics have many points of similarity [pT|,|8|,p^ ; 
in each case the low temperature phase is a gas of weakly interacting thermally populated 
stringy particles whose free energy is 0{1) in the relevant coupling constants {jj or gs). 
The phase transition is driven by a condensation of these particles, and is marked by an 
order parameter that winds around the thermal time circle (Polyakov loop or Euclidean 
winding string) . The high temperature phase has a free energy that scales like the inverse 
square of the relevant coupling. 

If we accept that stringy thermodynamics can be understood in terms of deconfine- 
ment, it follows that the relation p{E) ~ exp{E /Th) must break down in the energy range 
which is relevant for the high temperature deconfined phase. We can explain this break- 
down as follows. In the high temperature phase of asymptotically free gauge theories, the 
density of states is approximately that of a 3 -|- 1 dimensional free field theory with A^^ de- 
grees of freedom on a space of volume V, p{E) ~ exp(i?'^/^A^^/^y^/^). A crude estimate of 
the energy at the point of transition between the two phases may be obtained by equating 
these two formulas for p{E), yielding E ~ N'^VTfj. Note that this energy is proportional 
to A^^ ~ 1/fi's- Gravitational effects are large at these energies for every ^f^, no matter how 
small, explaining the failure of the free string estimate. 

As we have commented in §2.1, this transition may occur below the Hagedorn temperature. 
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So far we lack a quantitative demonstration that large gauge thermodynamics 
exhibits sharp stringy featureslll. The deconfinement transition in flat space occurs at 
strong coupling and has so far resisted efforts at quantitative analysis. In this paper we 
will make progress on this important problem by analyzing conflning gauge theories on 
compact spaces^ of size R rather than flat space, taking RAqcd ^ 1- The thermal 
behavior of such a theory may be computed reliably at all values of the temperature; we 
will demonstrate that it displays stringy features. 

2.4- Large N deconfinement on compact spaces 

As we have argued at the end of the previous subsection, there are several good 
motivations to study the thermodynamics of Yang-Mills theories on compact manifolds. 
The AdS / CFT correspondence forces us to address this problem. In a conflning theory the 
flnite size of the manifold cuts off the running of the coupling at the scale 1/-R, introducing 
a dimensionless coupling constant into the problem. 

In the case of a conformal theory like the A/" = 4 SYM theory the new scale permits 
nontrivial temperature dependence of thermodynamic quantities. In this subsection we 
will take a flrst look at the thermal behavior of gauge theories on such spaces. 

Clearly, the thermodynamical properties of an SU{N) gauge theory on a compact 
space of size R ^ 1/Aqcd closely resemble those of the theory in flat space. This re- 
semblance cannot be perfect at flnite A^, however, as the transition between conflned and 
deconflned behavior must be smooth rather than sharp for any system with a flnite number 
of degrees of freedom. In the N ^ oo limit, on the other hand, the number of degrees of 
freedom is inflnite even at flnite volume, and the deconflnement transition mimics its flat 
space counterpart more closely; in particular it remains sharp. 

Note that the Polyakov loop is no longer a good order parameter in the flnite volume 
theory because its expectation value vanishes for kinematical reasons (Gauss' law makes 
it impossible to put a single fundamental quark on a compact manifold, independently 
of the phase the theory is in). However, as we will explain in §5.7, it is possible to 
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A notable exception is the case of two dimensional gauge theories, which one can explicitly 



solve in the large A'' limit in many cases (see, for instance, |3C,3ll]). In this limit one obtains a 
Hagedorn-like spectrum with Regge trajectories. Moreover, there is also a stringy description of 
two dimensional QCD [ p2[ . 

The AdS/CFT correspondence, which establishes a duality between a gauge theory on a 
sphere and a string theory, motivates this study. 
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define related order parameters, either by looking at the norm of the Polyakov loop or 
by introducing an infinitesimal amount of fundamental matter. Alternately, one can use 
the second order parameter discussed in §2.2 above; the high and low temperature phases 
are sharply distinguished by the fact that the free energy scales as 0{N'^) and C(l), 
respectively, in these phases. It is still reasonable to write the effective Landau-Ginzburg 
free energy for this theory as the theory of a complex scalar field (V) on the compact space, 
and we will use a variant of this description in our analysis later in this paper. 

We have argued above that a large confining gauge theory, compactified on a space 
of size R ^ 1/ Aqcd, undergoes a phase transition at a temperature of order Aqcd- In this 
paper we will study the opposite limit RAqcd ^ 1- The gauge theory is weakly coupled 
in this limit (recall that the gauge coupling stops running at scale on a compact space) , 
and the thermodynamics may reliably be computed at all values of the temperature. We 
will find that the phase transition persists in this limit at a temperature of order 0{^). 
In the case of the d = 4 Af = 4 SYM theory on S^, the AdS/CFT correspondence has 



established [|ll| , p!3[] that, at strong coupling A, the theory undergoes a phase transition at 
a temperature of 0{^). We will establish that this phase transition continues to occur at 
weak coupling A at a temperature of the same order. 

The weak coupling nature of the phase transitions we describe in this paper (at 
RAqcd ^ 1 for confining theories, and at A ^ 1 for the A/" = 4 SYM theory) allows 
us to analyze them in detail. Quite remarkably we will find sharp signatures of stringy 
thermodynamics in these theories. In particular, we will demonstrate that the low temper- 
ature phase is always characterized by a Hagedorn growth in the density of states (see also 
0,^). We view this result as evidence for the existence of string duals for weakly coupled 
gauge theories on compact manifolds. 



3. Partition function for free Yang-Mills theory on a compact space 

We now proceed to directly analyze the thermodynamics of weakly coupled Yang- 
Mills theory on a compact space, beginning in this section (and the next two) with the 
free theory. For simplicity, we restrict here to spaces (such as S'^~^) for which all modes 
of the various fields are massive. In most of this section we will discuss the specific case of 
a U{N) gauge theory with fields in the adjoint representation, but it is easy to generalize 
our arguments to more general theories. In particular, our exact results for the partition 
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function in §3.3 will be completely general. As usual, the behavior of thermodynamic 
quantities in the canonical ensemble is governed by the partition function 

Z{(3) = Yl ^'''= J p{E)x^dE, (3.1) 

physical states 

where p{E) is the density of states and we define x = = e~^/^. As we discussed 
earlier, it is essential to keep in mind that the constraint from Gauss' law for a gauge 
theory on a compact space implies that the physical states over which we sum must be 
gauge-invariant . 

3.1. Two matrix harmonic oscillators give Hagedorn behavior for large N 

By expanding the fields into modes on the compact space (e.g. spherical harmonics on 
S^), our field theory may be viewed as a quantum mechanical system with infinitely many 
degrees of freedom (one for each mode). In the free theory, these degrees of freedom are 
decoupled harmonic oscillators, each in some representation of the gauge group. Before 
embarking on our detailed analysis of this theory, we present a simple demonstration that 
even restricting to two such oscillators in the adjoint representation of U{N) yields a theory 
that displays Hagedorn-like behavior in the large limitil. 

Thus, we consider two such modes, each with unit energy, created by operators A'^ 
and in the adjoint of U{N). Physical states are gauge-invariant and correspond to 
traces of products of these operators acting on the Fock space vacuum. Single-trace states 
of energy E are specified by a series of E A^'s and B^^s inside the trace. The number of 
such states satisfies 

2^/E < n{E) < 2^, (3.2) 

where the upper bound arises since for each of the E positions in the trace we may choose 
or S^, and the lower bound comes because cyclicity of the trace equates a given state 
to at most E — 1 others. In the large A^ theory all of these states are independent, so the 
density of independent single-trace states is Hagedorn-like : 

piE) = 1^ ~ E-e^-^, (3.3) 

where the Hagedorn temperature in this case is Th = 1/(3h = l/ln(2). 
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Hagedorn-like behavior in free large systems was also observed in 
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For finite A^, we cannot have Hagedorn behavior at aU energies, since the high temper- 
ature behavior must be that of a field theory (or quantum mechanics in the toy example 
described above) with finitely many fields. The departure from Hagedorn behavior at 
high energies comes about because the oscillators have only finitely many gauge-invariant 
degrees of freedom (of order A^^), so that only a number of traces of order A^^ give inde- 
pendent state^i. Thus, at energies of order A^^, trace relations will cut off the exponential 
growth of states, resulting in a cross-over to field theory (or quantum mechanics) behavior. 
In the limit of large A^, this cross-over becomes a sharp transition, as we will see explicitly 
in the next sections. 

3.2. Exact partition function for N = oo 

We now show that the exact partition function of our free gauge theory in the strict 
N = oo limit (where no trace relations exist) may be obtained by simple counting argu- 
ments, following 

We consider a model with m > 2 matrix-valued bosonic harmonic oscillators; the 
number m can be finite or infinite, as in the example of a compactified (i-dimensional gauge 
theory. Let the z*'^ oscillator have energy Ei. We encode the spectrum of oscillators in a 
"single-particle partition function" z{x) = - x^^ where the sum goes over all oscillatorsil. 
Examples of such single-particle partition functions for compactified 3-1-1 dimensional 
theories will be provided in §3.5 below. 

Following the discussion above, the partition function of single-trace states with k 
oscillators is given by 

Zfc = ^t^iL + (positive) (3.4) 
k 

where the factor ofl/k compensates for cyclicity in the trace, and the + (positive) acknowl- 
edges that this is an over-compensation in cases where the trace breaks up into repeated 
sequences of oscillators. Since in the large A^ limit a single-trace state can have any number 
of oscillators, we need to sum this result over k to find a partition function for single-trace 
states, 

oo 

ZsT = Zk = — ln(l — z(x)) + (positive). (3-5) 

fc=i 



For example, for a single N x N matrix, tic(A^^^) may be expressed as a combination of 
products of traces of lower powers of A. 

To be precise, this is the single-particle partition function for the U(l) theory, since in general 
the single oscillator states are non-physical. 
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Actually, it is not too difficult to promote (|3.5| ) to an exact expression which correctly 
accounts for repetitions inside the trace. This careful counting utilizes Polya theory and 
was already performed in [§|,|5|. For the convenience of the reader we reproduce it in 
appendix B. The exact result in the large N limit is 

^ST = -f;^ln(l-z(x^)), (3.6) 

where f{q) represents the number of positive integers which are not larger than q and are 
relatively prime to q. Note that the first term in ( p.6|) is precisely the term we explicitly 
wrote down in ( |3.5[ ). 

Using ( |3.6|) , we may now write the full partition function of our model, summing over 
states with arbitrarily many traces. The space of multi-trace states is simply the Fock 
space of single-trace states, so the full multi-trace partition function Z is easily obtained 
from ZsT (taking into account the fact that the single-trace states behave as identical 
bosonic particles), 

oo ^ 

HZ) = V -Zsrix^ = 



oo 



UO ^ U(J / \ 

n q 

n=l 0=1 

(3.7) 

= -E^(Ev'(^))Mi-^(x')) = 

fc=l q\k 
oo 

= -5^1n(l-^(x'^)), 
k=i 

where we have used the formula '^q^^'^iq) = k. It is easy to generalize this result to 
include also fermionic oscillators; denoting hy zb the single-particle partition function of 
the bosonic oscillators, and hj zp the single-particle partition function of the fermionic 
oscillators, we obtain 

oo 

ln(Z) = - ^\Yv{l ~ zb{x^) + {~1)^zf{x^)). (3.8) 
fc=i 

The inclusion of fermions does not change the qualitative behavior, so it is sufficient to 
assume that all oscillators are bosonic for the purposes of qualitative discussions. 
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Equation ( |3.8|) gives the exact partition function in the strict N = oo theory. It is not 
correct for the finite theory because we have assumed all traces to be independent .0 
However, as long as the thermodynamics is dominated by states for which the number 
of oscillators in a trace is small compared to A^^ (which will certainly be true for low 
enough temperatures), trace relations will be unimportant and (|3.8| ) should give a very 
good approximation to the correct behavior. 

In order to understand the dynamics of ( ^.7] ) we note some obvious properties of the 
"single-particle partition function" z{x), namely that z{0) = 0, z{x) is a monotonically 
increasing function of x, and that (since m > 2) z{l) > 1. It follows that the equation 
z{x) = 1 has a unique solution (which we will denote by x = xh) for < x < 1. If we 
denote the corresponding inverse temperature by (3h {e~^" = xh), then ln{Z{(3)) is well 
defined for (3 > (3h, but diverges like — ln(/3 — (3h) as (3 tends to f^H from above. For 
/3 < (high temperatures) Z{(3) in ( |3.7|) is ill defined. All these features are reproduced 
by the formula iH 

Z{(3) = j dEp{E)e-^^, (3.11) 

where p{E) ^ e^"^ at high energies. Consequently, our system of oscillators has Hagedorn- 
like thermodynamics with the inverse Hagedorn temperature I3h- Note that z{x) = 2x for 
the two oscillator model of the previous subsection, yielding xh = 1/2 and (3h = ln(2) in 
agreement with the discussion in that subsection. 

The fact that Z{/3) diverges for T > Th = l3]j^ indicates that Th is a limiting 
temperature if is set strictly equal to infinity. Of course, this cannot be correct at any 
finite A^, no matter how large. Indeed, a system of mN'^ harmonic oscillators (or even a 
compactified field theory for which m is infinite but there is a finite number of oscillators 



When the number of fields in a trace is of order A''^ trace identities relate various states 
that are assumed to be different in ( |3.7D . For example, any 2x2 Hermitian matrix M obeys 
tr(M3) = i(3tr(M)tr(M2) - iT{Mf). 
25 Let 

In(Z) = j dEp{E)e-^'' (3.9) 

and 

ZsT = j dEpsT{E)e-^^. (3.10) 



Then, p' {E) ~ pst{E) ~ e^"^ jE at high energies. Note that (U) and (|3.10|) match (U) and 
(|2.1| ), respectively, upon setting d = de// = 1 in those formulas. 
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below any fixed energy) can certainly be heated to any temperature. It must instead be the 
case that at high temperatures (T > Th), the condition for the validity of (|3.8|) (namely, 
that the number of fields contained in the typical state is much smaller than A^^) fails; 
we will see in §5 that this is indeed the case. Thus, in the high temperature phase, trace 
relations invalidate the analysis above, and traces no longer provide a useful basis for U{N) 
singlets in the high temperature phase. This is true at any finite and so it is also true 
in the large limit of principal interest to us. In the next subsection we therefore turn to 
a more generally applicable method of dealing with the U{N) singlet condition. 

3. 3. Exact partition function for free Yang-Mills theory 

In this section, we will derive an exact expression for the partition function of free 
Yang-Mills theory with arbitrary gauge group and matter content, on any compact space 
for which all modes of the various fields are massive. As we have discussed, the basic 
excitations of such a theory are single-particle states, corresponding to oscillators which 
arise by expanding the various fields into physical modes on the space. The general physical 
state is obtained by acting with arbitrary collections of these oscillators on the Fock space 
vacuum, with the constraint that the total state should be a singlet of the gauge group. 

To obtain an explicit expression for the partition function, suppose we have bosonic 
modes with energies Ei in representations Ri of the gauge group and fermionic modes with 
energies E'^ in representations R[ of the gauge group, ii Then, the partition function may 
be expressed as a sum over the occupation numbers of all modes, with a Boltzmann factor 
g-/3£' corresponding to the total energy, and a numerical factor counting the number of 
singlets (physically allowed states) in the corresponding product of representations. 



{# of singlets in sym'^^ (Ri) sym^^ {R2) ® ■ • ■ ® anti^'^ {R[) anti"^'^ {R'2) ® • • ■}. 



Note that in order to correctly account for particle statistics, we must symmetrize (an- 
tisymmetrize) representations corresponding to identical bosonic (fermionic) modes. An 
explicit expression for the group theory factor may be obtained by using the fact that the 
number of singlets in a product of representations is simply the integral over the group 

Prior to this subsection we have assumed that all modes are in the adjoint representation. 
We relax that assumption in this subsection. 



00 



00 



00 



00 




(3.12) 



20 



manifold of the product of characters for these representations, as reviewed in appendix 
A. Thus, we find 

z{x) = I [du] n { f: n I e ^^-''Knt^'.i^^) \ ■ (^-is) 

The sums in brackets are generating functions for the characters of symmetrized (antisym- 
metrized) products of arbitrarily many copies of a given representation Ri {R[). As we 
show in appendix A, these turn out to have simple expressions given by equation (A. 8). 
Using this equation, we find that 

Z{x)= y"mexpi [expi ^~^^^^\ ^^'^XR'XUn 

I i;m=l I I i;m=l 

(3.14) 

Finally, since the sums over i in the exponentials are simply sums over the single-particle 
modes, it is convenient to define bosonic and fermionic single-particle partition functions 
for each representation R, counting the single-particle states in this representation without 
the degeneracy coming from the dimension of the representation and without any gauge- 
invariance constraints, 



z^{x)= z^{x)= x^i (3.15) 

Ri=R R'.=R 

With these definitions, we may express our final result for the partition function as 

Z{x) = I [rft/]exp jj^ f; - [z^{xn + i-ir^'z^ixn] XR{Un\ ■ (3.16) 

V R m=l ^ J 

Thus, the partition function for an arbitrary free field theory (with no zero modes) on a 
compact space may be expressed as a single-matrix integral, with the action determined 
in a simple way by the single-particle modes of the theory. 

3.4- Generalization to arbitrary chemical potential 

Before considering specific gauge theories, we note that the result (|3.16| ) may be easily 
generalized to the case of a non-zero chemical potential. For a set of commuting charges Q"^ 
which commute with the Hamiltonian, we may introduce chemical potentials Ha = ^^(Qa) 
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simply by replacing x^* — > x^* Y[a^<^' ( |3.12| ), where Qf is the a^^ charge of the i^^ 
mode. The resulting partition function is 



Z(x, 



{Qa}) = [ [dU] exp I ^ i [z§ix^, {CD + {-ir+'z^ix^, {C})] XRiUn 



(3.17) 

where the generalized single-particle partition functions are now given by 

z«(x,{qJ)= Y1 ^""'U^al (3.18) 

Ri=R a 

3.5. U{N) gauge theories with adjoint matter on x 

For most of the remainder of this paper, we will focus on U{N) gauge theories on a 
sphere with adjoint matter. In this case, using XadjiU) = tr{U)tr{W) (where we denote by 
tr without a subscript the trace in the fundamental representation) , we obtain the unitary 
matrix model 

Z{x) = j [(it/]exp|^ -^(^^(a;"") + (-l)'"+^2F(x""))tr(C/"^)tr((C/t)'")| . (3.19) 

It is straightforward to work out explicit expressions for the single particle partition func- 
tions for various types of fields on a sphere of unit radius. In appendix B (see also 0,^) 
we find that for (3 + 1) -dimensional scalars, vectors, and chiral fermions on x IR, 

zs.{x) = ^^±1^, zy,{x) = ZF.ix) = (3.20) 

In a 3 + 1 dimensional free U{N) gauge theory with ns scalar fields, ny vector fields 
and np chiral fermions (all in the adjoint representation), we should use in ( |3.19| ) 
zb{x) = nszs^ix) + nvzv4{x), zf{x) = npZFAix). Note that all of these functions (and 
the single-particle partition functions in general) increase monotonically from 2; = at 
X = (zero temperature), and in all dimensions above 0-1-1 they diverge at x = 1 (infinite 
temper ature)0. 

For SU{N) gauge theories, the only difference (beyond the fact that U must have 
determinant one) is that every tr(t/)tr(t/^) in ( |3.19| ) should be replaced everywhere by 
{tr{U)tr{W) — 1), the character of the adjoint representation in SU{N). 



For a d dimensional field theory compactified on S'^ ^, z{x) diverges as 2J\f'^°^ /{I — x)'^ ^, 
where M'^"^ is the number of single-particle degrees of freedom. 
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4. Path integral derivation of the matrix integral and an order parameter 

In the previous section we computed the partition function of free Yang-Mills theory 
on a compact space A4 by counting gauge-invariant states. However, Z{T) may also be 
computed by the vacuum path integral of Euclidean Yang-Mills theoryll on x S"-]- where 
the circumference of the thermal circle is /3 = 1/T. In this section we will rederive 
the expression ( |3.19| ) for the partition function by computing this path integral for gauge 
theories on an of unit radius. 

Though it gives an identical final result for the partition function, the path integral 
derivation provides a physical interpretation for the unitary matrix U as the Wilson loop 
of the gauge field (averaged over the compact space) around the thermal circle. Con- 
sequently, j^tr{U) is precisely the Polyakov loop operator described in §2.2; a standard 
order parameter for deconfinement. We will see in the next section that -^tr{U) is also 
the natural order parameter for a large phase transition of the matrix model (|3.19| ) . 

A further advantage of the path integral derivation is that it generalizes easily to the 
calculation of the large partition function Z{T, A) at weak 't Hooft coupling A = Qym-^ 
when the theory is no longer free. In §6.2 we analyze the general structure of Z{T, A) at 
small A. An explicit calculation of the weak-coupling partition function, in specific gauge 
theories, is postponed to a future paper [|10| . 



4.1. Basic set-up 

We work in the gauge 



d^A' = 0, (4.1) 



where i = 1,2,3 runs over the sphere coordinates, and di are covariant derivatives. The 
choice ( [4.1|) fixes the gauge freedom only partially; it leaves spatially independent but time 
dependent gauge transformations unfixed. We fix this residual gauge-invariance with the 
condition 

dta{t) = 0, (4.2) 

where 

a = — [ Ao, (4.3) 

'^3 Js^ 



We thank Nima Arkani Hamed for emphasizing the utihty of this approach to us, and for 
several extremely enjoyable and productive conversations on related issues. 
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where ^3 is the volume of S^. 

The mode a wiU play a special role in what follows because it is the only zero mode 
(mode whose action vanishes at quadratic order) in the decomposition of Yang-Mills theory 
into Kaluza-Klein modes on x S^. Consequently, a fluctuations are strongly coupled 
at every value of A, including in the limit A — * 0. In particular, they cannot directly be 
integrated out in perturbation theory. 

In order to proceed with the perturbative evaluation of the partition function we will 
adopt a two step procedure. In the first step, discussed in this section, we integrate out 
all non-zero modeJH and generate an effective action Seffia) for a. Seff{(y) is non-trivial 
even at zero coupling and, as we will discuss in §6, it is further modified perturbatively in 
A. Once we have obtained S'e//(a), we must then proceed to perform the integral over a 
in order to obtain Z{T). This is the subject matter of §5. 

4-2. The integration measure 

In this subsection we will define Seff{ct) more carefully. Recall that the Yang-Mills 
free energy may be written as 



e-P^ = da DAAiAae-^^^^^'"^ (4.4) 



where Ai is the Fadeev Popov determinant conjugate to ( [4.1| ), A2 is the Fadeev Popov 
determinant conjugate to ( |4.2| ), and Sym is the Yang-Mills action. It is not difficult to 
explicitly evaluate A2 (see later in this subsection) and verify that it is independent of A. 
Consequently, (^4.4| ) may be rewritten as 



e 



-^^ = j daA2exp[-Seff{a)], (4.5) 



where 



exp[-5e//(«)] = j VAAiexp[-SYM{A,a)]. (4.6) 

In the rest of this subsection, we will explicitly evaluate A2 and so determine the effective 
measure of integration ( [4.5| ). 
It follows from (^) that 

A2 = det'(c>o-z[a,*]), (4.7) 



As a is a zero mode this is a standard Wilsonian procedure. Note the parallel with the 
discussions in §2.2 and in appendix C. 
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where the prime asserts that the determinant is over non-zero modes of Aq. Denoting by 
Xi (z = 1, ■ ■ ■ , n) the eigenvalues of a, and choosing a convenient basis of matrix functions 
whose time-dependence is given by exp(27rmt//3), the determinant is easily evaluated as 
the product 



n/O i,j 



27iin 

z [Xi - Xj) 



27vim \ T-r 2 f f3 (Xi - A,) 

III gjj^ I 



n^71 IIIL 1 T — r 



Notice that up to an overall constant, 

da A2 = [dU], (4.9) 

where 

da = l[dKl[{X,-Xjf (4.10) 

i i<j 

is the left-right invariant integration measure over Hermitian matrices a, and 

m=n^^^n^-'(^^^^) (4-11) 

i i<j 

is the left-right invariant integration measure in the integral over the unitary matrices 

U = e^^". (4.12) 

We will see in the next section at one loop order and argue generally in §6 that S^f / may 
be regarded as a function of U rather than a, so that ( [4.5| ) may be written as 

e-/^^ = j [dt/]e-^^//(^) (4.13) 

where Seff{U) is defined in (|4.6|). 

4-3. Evaluation of Seff at one-loop 



The path integral in (|4.5|) may be evaluated diagrammatically, generating an expansion 
of Seff{U) in powers of the gauge coupling. In this subsection we use this method to 
evaluate Seff{U) to lowest order in A ((9(A°)). This evaluation is rather simple as only 
one-loop graphs in the expansion of (^4.5| ) contribute to Sef / at this order. 

25 



We first consider the pure gauge theory. The Fadeev-Popov determinant for gauge- 
fixing is 

detdiD' = I VcVce-^^^^'^, (4.14) 



where denotes a gauge covariant derivative 

D,c = d,-i[A,,c] (4.15) 

and c and c are complex ghosts. Consequently, the action of the free gauge theory may be 
written as 

^-s,ff{u) _ j vA,VAoVcVcd{d^A') exp - ^ Tr (^^A,{DI + d^)A' + ^Aod'^Ao + cd^c^ 

(4.16) 

where 

DoX = doX -i[a,X]. (4.17) 
To proceed further we note that any vector field on the sphere may be decomposed as 

A, = d,^ + B, (4.18) 

where diB^ = 0. The integral over (f in ( [4.16|) is easily performed using the S function, 
yielding 1 / a/ det'{did^) where the derivatives act on scalar functions on and the prime 
denotes omission of the zero mode. The integral over Aq yields the identical factor (the 
zero mode is a which is not integrated over). The integral over the ghosts, on the other 
hand, evaluates to det'{did^). These three factors cancel nicely, so that ([4.16|) simplifies to 



JvB^exp -^jTv(^B,{Dl+d^)B' 



(4.19) 



Thus, 



5e// = ^ln(det(-D2-a2)), (4.20) 



where the operator acts on the space of divergenceless vector functions on the sphere, i.e. 
the space of vector functions spanned by the vector spherical harmonics. Consequently 



Seff = lYl ^(^) lndet(-I^2 + A'), (4.21) 



2 A 



where are the eigenvalues of the Laplacian acting on vector spherical harmonics on 
the compact manifold and n(A) is the degeneracy of each eigenvalue. When the compact 
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space is an of unit radius we have A = h + 1 (for integer h > 0) and n(A) = 2h{h + 2) 
(see the counting of vector spherical harmonics in appendix B). 

The determinant of (— -Dq + easily evaluated by passing to Fourier space in the 

time direction, yielding the infinite product 



det 



U{N) 



n 



n=— oo 



— + — " + ° +^ 



(4.22) 



The infinite product of matrices may be rewritten as 

-2 r 



n 



27rm 



n -2 



n 







27rm 



n 

-2 



/3 



2Tvm 



n=l 



/3^(« + zA)- 

4yj-2^2 



/3^(a-iA)- 

4yj-2^2 



2 



n 



/9 



27rm 



2 



sm 



/3(a + iA) 



sm 



f3{a - iA) 



[cosh(/3A) - cos(/3a)] 



where 



A/" 



v 



1 



n 



27im 



(4.23) 



(4.24) 



The divergent factor A/" is a constant independent of both A and a, and we set it to 
unity to reproduce the free energy of the harmonic oscillator. Thus 



ln(det(-L'2 ^ ^2^^ = Tr (/3A + ln(l - g-^^+^/^a^, ^ ^^^^ _ g-/3A-i/3a^)^ 



(4.25) 



where the trace is over A^^ x A^^ dimensional matrices and a acts in the adjoint representa- 
tion. Expanding the logarithms in a power series, summing over A, and passing from the 
adjoint to the fundamental (Trad, (e*"^^") tr(t/'")tr(t/'~")) and using ( [4.201 ) we obtain 



1 ^ / n\ 

Seff = -3N^Y.^<^) - E ^^tr(f/-)tr(f/--). 

ft 



(4.26) 



n=l 
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For the case of an of unit radius the first term (appropriately regularized) is equal to 
j^/3A^^, where is the Casimir energy J2o^ h{h + l){h + 2) for a vector field on the unit 
sphere.i^ 

It is a simple matter to generalize this calculation to include the contributions from 
free conformally coupled scalar and spinor fields. When the compact manifold is S^, the 
contribution to Seff from a single additional scalar field is 

5S,ff = ^ ln(det(-52 -d^ + 1)), (4.29) 

where the operator acts on scalar fields on x and the constant piece of the operator 
is a consequence of the TZf'^ term in the Lagrangian for conformally coupled scalars. As 
above, this determinant is easily evaluated to yield 

^^ff = 240^^' ~ S ^^tr(f/-)tr([/-), (4.30) 

n=l 



See equation (64) of |34]. This equation, and the other Casimir energies in this subsection, 
may be justified as follows. On a sphere of radius R we wish to compute X]m-i('^^ ~ 1)('^/^)) 
where E = m/R is the energy of the mode in question. We regulate this sum by multiplying 
the summand with a cutoff function f(E/A), where / is a smooth function such that /(O) = 1, 
/'(O) = and /(oo) = 0. We now evaluate the regulated sum using the Euler-MacLaurin formula 

-F(0) + F(l) + F{2) ... = j F{x)dx + C(-l)i^'(0) + C(-3)^P + • • ■ (4.27) 

with F(m) = (m^ - l)(m/i?)/(-^), and we find 

K'R^ x'fix) - K'rJ^ xf{x) + + + (4.28) 

In order that the vacuum energy of the theory be zero in flat space (a necessary condition, for 
instance, for conformal invariance), the two divergent terms above must be cancelled by counter- 
terms; indeed, counterterms of the form aA"* J ^ and bA^ J (with suitable values for a and 
h) achieve this cancellation. The remaining finite piece of the energy is 11/120-R. Finally, note 
that we have been careful to choose a 'general coordinate invariant' regulator; one that cuts off 
modes in a manner that depends only on their proper energy. It would be incorrect, for instance, 
to regulate the summand above with the smoothing function /((m — 1)/Ai?); in this regulation 
scheme theories with the same Lagrangian and same A but different R are not identical in the 



UV. See chapter 6 of [35| for an interesting and extensive discussion of related issues. 
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where = | X]^o('^ + 1)^- A Weyl spinor field on 5"^ contributes 



S,/; = -ln(detMo7''-?)) = 

960 n 



where ^ = - Er=i ^(^ + + l)- 

Thus, we have rederived the expression for the free partition function ( |3.19| ), but this 
time with a physical interpretation for the matrix U as the holonomy around the time 
circle, and for -^tr(t/) as a Wilson loop around the time circle, or a Polyakov loop (in the 
fundamental representation) . 



5. Solution of the free Yang-Mills matrix model 

In this section we proceed to directly analyze the unitary matrix model ( |3.19|) for the 
exact partition function of free U{N) Yang-Mills theory in order to extract the thermody- 
namic behavior. 

We begin by recalling that any unitary matrix model with gauge-invariant action and 
measure may be rewritten entirely in terms of the eigenvalues of the unitary matrix, which 
must lie on the unit circle. Denoting these eigenvalues by {e*"^} (with — tt < < tt), we 
may rewrite the partition function (|3.19D in terms of the eigenvalues by the replacements 



/ m [da,] n sin^ (^^) ; tr(C/-) ^ Yl ^^^"^ ' (^-l) 

i ^ i<j ^ ^ j 

With only adjoint matter, we find from ( |3.19| ) that the effective theory for the eigenvalues 
is governed entirely by a pairwise potential, 

Z{x)= /"[rfa,]e"^»^^^^"'""^^ (5.2) 



where 

oo ^ 

V{e) = - In I sin(6'/2)| -J^- [zb{x'') + {-ly+hpix'')] cos(n6') 



n 

(5.3) 



ln(2) -f V -(1 - 2b(x") - (-ir+^Fixn) cos(n6) 



n 

n=l 



We take the spinor field to be anti-periodic on the in order to compute the trace of e 
rather than the trace of ( — l)^e^^^. 
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In the first line, the first term coming from the measure is a temperature-independent 
repulsive potential, while the remaining terms provide an attractive potential which in- 
creases from zero to infinite strength as the temperature increases from zero to infinity.il 
This suggests that at low temperatures, the minimum action configurations will corre- 
spond to eigenvalues distributed evenly around the circle, while at high temperatures, the 
eigenvalues will tend to bunch up. 

For finite N, the partition function receives contributions from all configurations of 
eigenvalues and depends smoothly on the temperature. On the other hand, it is well 
known that in the large limit, the matrix model partition function is dominated by the 
minimum action configurations, with the exact leading and subleading terms in the 1/N 
expansion of the free energy given, respectively, by the minimum of the action and by the 
Gaussian integral about the minimum action configuration. When this minimum action 
configuration (or its derivatives) changes abruptly as a function of temperature, we may 
have a phase transition in the large limit (see, for example, 0). We will now see that 
exactly this behavior occurs and leads to a phase transition for our model. 

We begin in §5.1 by analyzing the low temperature phase. In §5.2 we analyze the be- 
havior near the phase transition temperature, and in §5.3 we analyze the high temperature 
limit. In §5.4 we write down the solution of the matrix model in the high temperature 
phase, and in §5.5 we provide a perturbative expansion of this solution around the phase 
transition temperature. In §5.6 we summarize the thermodynamical behavior of the the- 
ory. In §5.7 we discuss how the Polyakov loop may be used as an order parameter at finite 
volume, and we end in §5.8 with an application of our results to some interesting 3 + 1 
dimensional gauge theories. 

To see that the second term in the potential is always attractive, note that we may rewrite 

it as 

V2ie) = ^ j dE[pB{E)\n{l-2x'^ cos{e) + x^^)- pF{E)\n{l + 2x'^ cos{e) + x^^)] , (5.4) 

where pB and pp give the single-particle density of states for bosonic and fermionic modes . This 
potential is always attractive since for any value of E and x the integrand is decreasing in the 
interval 6 G [— tt, 0] and increasing in the interval 6 G [0,7r]. 
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5.1. Low temperature behavior 

Consider first the theory at low temperatures, where the effective potential between 
the eigenvalues is dominated by the repulsive term. In this case, we expect that the 
pairwise repulsion drives the eigenvalues to spread uniformly around the circle. In fact, 
since there is no difference between displacing any individual eigenvalue in the uniform 
distribution to the left and displacing it to the right, the uniform eigenvalue distribution 
will always be a stationary point of the action for a pairwise potential. To see when this 
stationary point is a minimum, it is convenient to introduce an eigenvalue distribution p{9) 
proportional to the density of eigenvalues e*^ of U at the point 9. Note that p must be 
everywhere non-negative, and we may choose its normalization so that 



dOpie) = 1 . (5.5) 
With this definition, the effective action for the eigenvalues becomes 
S[p{9)] = N^ [ dOi [ d92p{ei)p{e2)V{9i - 62) = 



27V 



(5.6) 



n=l 



where in the second line, we have defined p^ ^ f d6p{6) cos{n6) and = f d6V{6) cos{n9) 
to be the Fourier modes on the circle of p and V, respectively, and we assume without 
loss of generality that the eigenvalue distribution is symmetric around ^ = 0. From the 
latter expression, it is clear that the uniform distribution {pn>i = 0) will be an absolute 
minimum of the potential as long as all are positive. From (|5.3|) we see that 



K = —(1 - zb{x^) - {-ir+'zp{xn), (5.7) 
n 

so the uniform distribution is an absolute minimum if and only if 

ZBix"") + (-l)"+izi.(x") < 1 (5.8) 

for all n. But, since the single-particle partition functions are monotonically increasing, 
and < X < 1, the n = 1 condition is always strongest, and the uniform distribution will 
be stable for temperatures T < Th = —1/ ln{xH) where xh is the solution to 

z(xh) = ZBixn) + zf{xh) = I ; (5.9) 
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note that this is precisely the Hagedorn temperature we discussed in §3.2. As long as we 
have more than a single quantum- mechanical mode, equation (5^) always has a unique 
solution with < x < 1, so the uniform distribution becomes an unstable extremum 
beyond some finite temperature Th determined by the single-particle partition function. 

For T < we may now evaluate the free energy at leading and sub leading order in 
A^. Since tr(t/"^) vanishes for any n > 1 for the uniform distribution, the classical value for 
the action (|3.19| ) (and thus the leading 0{N'^) contribution to the free energy) vanishes. 
The first non-zero contribution to the free energy thus arises from the Gaussian integral 
around this configuration. From the quadratic action (|5.6| ),( [577D , we see that the leading 
contribution to the partition function is ii 

oo ^ 

n=l \ / \ / V / 

where the normalization has been arbitrarily fixed by choosing the vacuum free energy to 
vanish, Z{x = 0) = 1 (more physically, we could choose Z{0) to account for the Casimir 
energy of the vacuum, see the previous subsection). This is precisely the same as the result 
( ^.8| ) that we found above by counting the states with E -C A^; we see that this is indeed 
the leading large A behavior as long as T < Th- 

This partition function diverges at T = Th, as should be expected, since this is the 
point where the quadratic action ( |5.6| ) develops an unstable direction. The corresponding 
divergence in the free energy F = — Tln(Z) goes like 

F^ThHTh-T) (5.11) 

as T approaches Th from below. By the general discussion in §2, we may conclude that 
this divergence is associated with a Hagedorn density of states, 

p{E) oc EV^-^ , (5.12) 

as we also found by counting states in §3. Thus, just as for perturbative string theory, we 
have a Hagedorn divergence which is associated with an unstable mode in the Euclidean 
path integral; note that the power of the energy is also the same as in string theory at finite 
volume. In this case, the "tachyon" causing the divergence is the lowest "momentum" fiuc- 
tuation mode pi of the eigenvalues. In the next subsection, we will be able to see explicitly 
what the endpoint of this tachyon condensation brings as we raise the temperature beyond 
the Hagedorn temperature. 



The 'Jacobian' for the 'change of variables' from Ai to pn appears to be irrelevant at every 
order in the 1/N expansion, see ||3^ . 
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5.2. Behavior near the transition 

At the temperature T = Th, Vi vanishes and the mode pi (the lowest momentum 
fluctuation in the eigenvalue distribution) becomes massless. Since the action is quadratic, 
this corresponds to an exactly flat direction in the potential, and the corresponding family 
of minimum action conflgurations are given (up to overall translations of 6) by 

p{e) = ^{l+t cos(0)), < t < 1, (5.13) 

ZTT 

which interpolate between the uniform distribution and a sinusoid vanishing at a single 
point. For T > Th, the quadratic form S has negative coefficients, so all minimum 
action conflgurations must lie at the boundary of conflguration space, at a point where a 
hyperboloid of constant S lies tangent to this boundary. Since this boundary is provided by 
the positivity condition p{9) > 0, these minimum action eigenvalue distributions necessarily 
vanish on a subset of the circle. 

In the limit of small, positive AT = T — Th, the action contour 5" = is a cone with 
opening angle going to zero. The contours of smaller S are hyperboloids inside this cone, 
so the minimum action conflguration lies on the boundary of the conflguration space inside 
the cone. Finally, it is easy to see that the region of allowed p's is convex (if points {pn} 
and {p'n} are in the region then all points {tpn + (1 — t)p'^} are in the region for < t < 1). 
Thus, the boundary region interior to the cone is a simply-connected neighborhood of the 
t = 1 conflguration in (|5.13D whose size goes to zero as T — >^ Th- We conclude that the 
minimizing conflguration changes continuously from the t = 1 conflguration in (|5.13|) for 
T in some interval above Th- 

The leading behavior of the free energy in this vicinity of Th is then given by evaluating 
the Vi term in the action ( ^.6|) on the t = 1 conflguration in ( ^.13|) , since the effects of 
changing the eigenvalue distribution come in at higher orders. 

Thus, the leading result for the minimum action is 

St>tJ^) = ^(T - Th) VHTh) + ... (5.14) 
from which we may obtain the free energy as F — TS"^"^^. The leading behavior near the 
transition is therefore given by 

so we have a flrst order phase transition at the Hagedorn temperature. Note also that 
the behavior is characteristic of a deconflnement transition, since the free energy is 0{1) 
below the transition and 0{N'^) above the transition. This will be verifled in §5.7 below 
when we discuss the behavior of a second order parameter for conflnement, related to the 
Polyakov loop. 
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5.3. High temperature behavior 

For general values of T > Th, we will not be able to write down an exact expression 
for the free energy, though we describe in the next subsection a formal solution that may 
be used to evaluate the free energy to arbitrary accuracy. However, the theory simplifies 
again in the large temperature limit, where the potential becomes strongly attractiveil. In 
this case, the action is minimized by a tightly clustered configuration of eigenvalues that 
approaches a delta function in the limit of infinite temperature. Thus, we should have 
= 1 to leading order, and 

5=— 5^K(T»1). (5.16) 

n=l 

The high temperature limit of the single-particle partition functions depends only on the 
space-time dimension d and the number of physical polarizations J\f'^°f of the fields in the 
theory, 

z^{x) ^ 2Mt"^T'^-^ + 0(T^-2) . (5.17) 



Using (|5.7|) , we find that the high temperature limit of the free energy is 



F{x ^ 1) = -2N^T'^C{d) 



<^+(l--^X°^ 



+ 0{T'^-^). (5.18) 



The limiting free energy density here coincides with that of the flat space theory, which 
should be expected since in the free field analysis we are dealing with a scale-invariant 
theory, so taking the dimensionless temperature TR to infinity is equivalent to taking 
the limit of large volume at fixed temperature (see [ P5| , |5D|J^ among others for the high 
temperature expansion of the free energy for the U{1) theory). 

5.4. Exact solution for T > Th 

Having understood the qualitative features of the T > Th behavior, we now describe 
an exact solution for the matrix model at arbitrary temperatures above the transition. Just 
beyond the Hagedorn temperature, we have seen that the eigenvalue distribution becomes 
sinusoidal, vanishing at ^ = tt as the lowest momentum fluctuation mode condenses. As the 
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The high temperature regime and its relation to string theory was discussed for uncompact- 



ified theories in [37| in a very similar language to the one we are using here, and this discussion 



was applied to two dimensional QCD in |31]. 
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temperature is increased further, the attractive term in the pairwise potential continues 
to increase in strength, so the eigenvalues will become increasingly bunched together, 
occupying only a finite interval / = [— ^o, ^o] on the circle (we arbitrarily choose the 
middle of this interval to be at ^ = for convenience). 

The precise distribution may be determined by the condition that a single additional 
eigenvalue added in the interval / experiences no net force from the other eigenvalues, 

/ V'{a~e)p{e)de = 0, ael. (5.19) 



Substituting in our potential (|5.3|), we obtain 



'0 



cot 



p{9)d9 = 2 ^ ttnPn sin(na), (5.20) 



\ / n=l 



where on the right hand side, a„ = zb{x^) + {—1)^~^^zf{x^), and we have used the fact that 
the eigenvalue distribution should be symmetric about ^ = and therefore be orthogonal 
to sin(n^) for all n. As before, pn are the moments of the eigenvalue distribution. 



f p{9)cos{n9)d9. (5.21) 

J —n 



To solve for p{9), it is convenient treat the moments p„ as independent variables, and 
self-consistently solve for p{9) and pn together using the two equations ( p.20|) and ( |5.21|) . 
Fortunately, the equations ( |5.20| ) provide the equilibrium conditions for a matrix model 
with action 

oo 

anPr. 

n 



5 = AT J]^I^(tr(C/")+tr(C/t")), (5.22) 



n=l 



that has been solved exactly in [^1[. The solution for p{9) i 



p{9) = i^sin^ (I) -sin^ cos((n - ^)^) (5.23) 
for —6*0 < 9 < 9q and p{9) = otherwise, where 

oo 

iPiicosi9o)) (5.24) 



1=0 
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In the language of the solution we are looking for is of type Ai. 
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and the angle 6*0 that bounds the eigenvalue distribution is determined by the condition 



gi=go + 2. (5.25) 
In the formula (|5.24| ) for Qn, the Pi are Legendre polynomials, defined by 

oo 

^Pi{x)z^ = {1-2XZ + z'^)-^. (5.26) 



1=0 



We may now eliminate p{9) to obtain a linear system of equations for the pn from ( |5.21j ) 
and ( ^.25|) . To describe these, we define a matrix R and a vector A whose elements are 
polynomials in = sin^ (6*0/2), 



R^i ^ ai 5^(S"^+'=-i(s2) + Sl— ^+il(s2))p^_,(i _ 2s' 
k=i 

Am = am{Pm-i{l - 2s^) - Pm{l - 2s^)), 

where 



(5.27) 



B—-{s') = -J^ de^sm' ^-fj - sin^ ^- j cos((n - -)e) (5.28) 
are polynomials defined by the generating function 

^S"+5(a;)2" = —{y/{l-z)'^ + 4zx + z-l). (5.29) 

n=0 

The determining equations ( [5.21| ) and (|5.25| ) for the eigenvalue moments and the angle 9q 
are then simply 

Rp = p; A-p=l. (5.30) 

Thus, the vector of moments of the eigenvalue distribution must be an eigenvector of the 
matrix R with eigenvalue 1, normalized so that its dot product with A is 1 (ensuring that 
the eigenvalue distribution integrates to 1). The condition 

det(i?-l)=0 (5.31) 

that the matrix R has such an eigenvector determines the angle 9q. Note that while this 
equation is a complicated function of 9o, it is linear in each of the coefficients a^, so it is 
convenient to parameterize the solution in terms of (sin(^o/2), {an>i})7 solving for ai in 
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terms of these variables using (|5.31 ). With ( |5.31|) satisfied, the exphcit solution for the 



unit eigenvector may be given as 

p = M-^ei, (5.32) 

where M is a matrix obtained by replacing the first row of the singular matrix (1 — i?) 
with the vector A, and ei = (1, 0, 0, . . .). 

This solution is rather formal for the general case when R is an infinite matrix. How- 
ever, in the case when an>k vanishes for some k, the matrix M takes the form 

^=(^l' l)' (^-^^^ 

so the explicit solution may be obtained by inverting a finite matrix. 

Note that once we evaluate the moments Pn<kj the full eigenvalue distribution is given 
by (|5.23|) since Qn>k = in this case. Similarly, only the upper k x k submatrix of R 



contributes to the condition ( ^.31|) determining ^o- 

In our model, the coefficients die off with n (exponentially for a; <^ 1, and like a 
power law as x — > 1) so we may obtain an arbitrarily good approximation to the exact 
eigenvalue distribution by truncating such that an>k = for sufficiently large k. It turns 
out that even restricting to /c = 1 retains the same qualitative behavior of the model.ll In 
this case, the angle 9o is determined by 

= det(i?ixi - 1) = ai(2s^ - s"^) - 1, (5.35) 

while pi is determined by 

pi = M-^^ ■ 1 = {2ais^)-\ (5.36) 

The full distribution is then determined from pi using ( [5.23| ) (noting that Qn>i = 0). 
Parameterizing the solution by s = sin(^o/2), the final result is 




P{0) = — ^^^T^Wsin' ( ^ ) - sin' ( o ) ^os ( - ) , (5.37) 
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In this case, the auxiliary theory ( ^.22| ) is exactly the model studied in Q. The truncation 
to A; = 1 was also considered in H. 
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with 6*0 determined in terms of ai by (|5.35|) , or explicitly 



2 / ^0 , . 
sm I y 1 = 1 - 



ai(T) 



(5.38) 



As T increases from Th to infinity, ai increases from 1 to infinity, so decreases from tt 
to 0, such that the eigenvalue distribution (|5.37| ) eventually approaches a delta function. 
Using (|5.37|) , we may evaluate the free energy in the approximation a^yi = and find 
that 



1 

'iV2 



ln(Z, 



t rune J 



1 1 w 2x 1 

+ - In(s^) 



(5.39) 



Near the transition, we find 



+ 0((T-Th)'), (5.40) 



T=Tn 



which gives the same leading behavior ( [5.15| ), characteristic of a first order transition, as 
we found in the full model. 



5.5. Perturbative expansion slightly above the Hagedorn temperature 

Using the general solution above, it is also possible to determine explicitly the exact 
behavior near the transition as a perturbation expansion in (T — Th). To determine 9o as 
a function of temperature, it is convenient to define A = ai — 1 and e = cos^ (6*0/2), and 
to expand R = Rq + SR where 

i^o = diag(ai, a2, as, . . .), (5.41) 

so that 5R will be small for small (T — Th). We may then solve the condition ( |5.31J ) 
perturbatively for A as a function of e by writing 

= det{l-{l-Ro)~^SR) (5.42) 

and expanding the determinant out in terms of traces. This expansion in powers of 5R 
may further be expanded as a power series in e and A, and finally solved perturbatively 
for A as a power series in e. Using the fact that 

SRmi = {-l)^+'+H^maie^ + 0{e^) , (5.43) 
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and that the lowest power of A appearing in ( |5.42| ) is 1/A (from (1 — Rq)^^), it follows 
that A may be obtained to order e"^^^^ by keeping terms to order SR^ in ( p.42|) . Explicitly, 
we find that the leading terms are 



A = e' 



n Gr 



(5.44) 



n=2 



When a„>i = 0, this correctly reduces to the leading orders of the truncated model result 

(S) 

n>l 



To find perturbative expressions for the moments of the eigenvalue distribution near 
the transition, we use the general result (|5.32|) , expanding M = Mq + 5M with 



Mn 



1 - a2 

1 - as 

V -J 



/4 



Ai Ai(l-a2) Ai(l-a3) 
1 

1 — a2 



1 — as 



SO that the explicit perturbative solution is 



\ 



(5.46) 



p = Mo\l + SMM^^y^ei = (Mq-^ - Mq^SMMq^ + . . .)ei . 



(5.47) 



Note that 



ml 



-SRml 




m > 1, 
m = 1, 



(5.48) 



so from ( |5.43|) we see that SM starts at order e^. Thus, we may obtain the eigenvalue 
moments to order e^fc-i |-,y keeping the first k terms here. Using the first two terms shown 
and the expansions of ai and in powers of e, we find 



Pi = (Mo-i)n + J2(^o')inSRniiM^')ii + 0{e' 



n=2 

1 e 



Lee \ nan 

= - + - + —> — + 0( 

2 2 2^1- a„ ^ 



n=2 



(5.49) 



p^>i = {M-XnSRniiM-^)n + Oie^) 



n 



'-ir—— + 0{e'] 
1 - a„ 
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The result for Pn>i seems at first somewhat surprising since it appears that the moments 
increase without bound as n is increased. However, it turns out that the leading terms 
( p.43| ) in the expansion of SRni in e provide a good approximation only for e ^ 
To obtain an approximation valid for all n at some fixed value of e, we should use the 
complete expression 6Rni = ai(S"+2 (1 - e) + ^""2 (1 - e)), which does fall off for large 
71. Alternately, we note that the perturbative expansion of the difference between the 
moments of the exact distribution and those of the truncated k = 1 model is well behaved 
for large n, 

6pi = Pi- p,{a^>, = 0) = V E TZ^ + ^(^')' 

Z _ 1 CLm. 



(5.50) 



5pn = Pn>l - Pn>l(an>l = 0) = -("1)":^ — + 0{t 

2 1 - a„ 



due to the exponential decay of a„ for large n. 

Using the expressions ( ^.49|) for p^, we may finally evaluate the free energy from (| 
and ( p77|) . To avoid subtleties associated with the apparent growth in the moments ( |5.49| ) 
for large n, it is again convenient to express the result as a correction to the truncated 
model free energy. The final result, taking into account all corrections to the truncated 



model, is 



1 °° 1 

Jp ln(Z) = X] -(1 - an){pn{an>l = 0) + 5pnf = 

1 , X e'^ v-^ niin? — l)a„ =, 

-J^HZtrun.) + \_J + ) = (5-51) 



n=2 



where Ztmnc is the partition function ( ^.39|) of the k = 1 truncation computed earlier in 
this section. If desired, the free energy may be rewritten in terms of A using (|5.44|) , or 
explicitly in terms of temperature using 

A = ai(T) - ai(TH) = a[{TH){T - Th) + \'i{TH){T - Th? + • ■ ■ . (5.54) 



"^^ While the perturbation expansion about the k = \ truncated model in the second line of 
( 5.51| ) is perfectly well behaved, the naive evaluation of ln(Z) using ( 5.49D yields the divergent 
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lim (F/NO 

N 



Figure 1: Free energy as a function of temperature in free Yang-Mills 
theories. 

5. 6. Summary of thermodynamic behavior 

We are now in a position to present a reasonably detailed picture of the behavior of 
the free energy of the free gauge theory at all temperatures, depicted in figure 1. 



expression 



1 ^1 
■— ln(Z) = ^-(l-a„)p: 



n=2 



n 



2 + 4 



n=2 



n(n 



l)ar. 



(5.52) 



The divergence in the third term of ( |5.52D is unphysical; it is a consequence of the (apparent) 
hnear growth in the moments observed in ( ^.49| ). As we noted earher, this hnear growth is actually 
cut off beyond n ~ where it is replaced by a decaying behavior. Thus, the divergent sum in 

( 5.52| ) should be cut off at n = C / ^/t, resulting in a contribution that scales like a finite 

contribution to a lower order of perturbation theory. Indeed, it is easily verified that the (finite) 



term in (|]5|) disa grees with the exact term in the expansion of the solution (|5.39| ) 



-—\n{Ztr. 



3e^ 



0(e^ 



(5.53) 



as our analysis above would suggest. We conclude that the divergent sums at higher orders in 
the naive perturbation theory conspire to give finite contributions + e'^/S + • • •. A systematic 
understanding of this naive perturbation theory would be an intricate task; one that may, however, 
be avoided by choosing the correct starting point for the perturbation expansion, as we have done 
above. The smarter perturbation expansion ( 5.5l| ) is entirely divergence free; this follows since 
the a„ decay exponentially with n, whenever x <1. 
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We find that the uniform eigenvalue distribution provides an absolute minimum of 
the matrix model potential for all temperatures below some Th determined by the single- 
particle partition function as 

z{Th) = 1 . (5.55) 

The free energy in this regime, given by (|5.1C1| ), is 0{1) and decreases monotonically to a 
logarithmic divergence ( |5.11| ) at T = Th, characteristic of a Hagedorn density of states 

(iH). 

As T increases past Th, the uniform eigenvalue distribution develops an unstable 
mode which condenses to give a pure sinusoid distribution ( ^.131) for T just above Th, and 
a distribution with a gap for T > Th- For small {T — Th), the 0{N'^) free energy decreases 
linearly in T — Th , so the first derivative of the free energy is discontinuous and we have a 
first order phase transition. The free energy as a function of temperature may be written 
as a perturbation expansion in T — Th with the leading terms given in ( |5.51|) . 

In the high temperature limit, the pairwise potential between eigenvalues becomes 
strongly attractive, and the eigenvalue distribution approaches a delta function. In this 
regime, the free energy asymptotes to ( |5.18|) which corresponds to the expected flat-space 
free-energy density. Using the formal exact solution of the previous section, it should be 
possible to work out the high temperature behavior as a perturbation expansion in 1/T, 
as we did explicitly for the behavior near the transition. 

5. 7. The Polyakov loop as an order parameter at finite volume 

In the theory we analyzed above, the Polyakov loop (-^tr([/)) vanishes in both the 
high temperature and the low temperature phases. However, the reason for this vanishing 
is rather different in the two phases. The low temperature phase is governed by a single 
saddle point (uniform distribution of eigenvalues); the Polyakov loop vanishes on this saddle 
point. On the other hand the high temperature phase is governed by a one parameter set 
of saddle points (related to each other by the U{1) rotation U e'^'^U , which acts on the 
eigenvalue distribution by shifting 9 periodically). The value of (-^tr([/)) is non-zero when 
evaluated on each saddle point individually but vanishes upon integrating over a (averaging 
over all saddle point^l) , since the phase varies over a complete circleli. To obtain an order 
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Witten 1 11] has described a rather similar phenomenon in the bulk dual to the strongly 



coupled A/" = 4 Yang Mills. 

"^^ For SU {N) theories this is not an integral but a sum over A'' discrete saddle points, but it 
has the same effect. 
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parameter which correctly distinguishes the two phases, it is therefore sufficient to use the 
squared magnitude of the Polyakov loop, (iPp) = liniN^ooi^l^^iU)]'^) . 

This order parameter may be given a nice physical interpretation by the following 
alternative definition. Define {jjtr{U))e as the expectation value of {jjtr{U)) upon per- 
turbing the effective action in ( |3.19| ) by the infinitesimal term Ne{tic{U) + tr([/^)). Note 
that this perturbation breaks the U{1) symmetry (or the '^n symmetry for SU{N)); in 
the presence of this perturbation, the high temperature phase is also governed by a single 
saddle pointS. Consequently, the order parameter 

V^^f = lim lim (— trft/)), (5.56) 

e^O N^oo N 

is zero in the low temperature phase but non-zero in the high temperature phase0. It is 
easy to show that in the large limit, {\V\'^) = \V^^^\'^^ so the two definitions of the order 
parameter are equivalent. 

The construction presented in the previous paragraph may be understood physically 
as follows. As discussed in §2.4, (tr(t/)) vanishes in both phases for purely kinematical 
reasons; Gauss' law does not allow the introduction of a quark on a compact space. How- 
ever, the perturbation described in the previous paragraph effectively adds a condensate 
(or plasma) of classical quarks and anti-quarks to the theory; in the presence of such a 
condensate the constraint from Gauss' law is circumvented (recall that flux lines decay in 
a Higgs phase) . In the high temperature phase the additional classical quark has finite free 
energy even in the limit that the classical Higgsing is taken to zero. On the other hand, 
in the low temperature phase the Polyakov loop vanishes as this Higgsing is taken to zero. 

The analysis above shows that V^^^ = 1/2 immediately after the phase transition, 
and it goes up to V^^^ 1 at high temperatures. In fact, V^^^ is exactly the variable pi 
discussed in §5.4, so the exact behavior of the order parameter is given formally by the 
first component of p in ( |5.32| ), while the behavior near the transition is given by ( 5.47| ), or 
explicitly to leading orders by ( |5.49| ). 



Of course at any finite N there is some finite eo{N) that is so small that when e < €o{N) 
the path integral receives significant contributions from the entire manifold of 'almost' saddles. 
However to{N) — > as — > oo. Consequently, if we take the N ^ oo limit first, only a single 
saddle point contributes to the path integral for any e no matter how small. 

Note that this construction is rather similar to the definition of the order parameter of a fer- 
romagnet in 3 dimensions. Naively, the magnetization of the ferromagnet vanishes due to SO {3) 
symmetry in both the high and low temperature phases. However if we define |M|<: to be the mag- 
netization in the presence of a small magnetic field in the z direction, then lime^o limy^oo \M\e 
(where V is the volume) is the standard order parameter for the ferromagnet. 
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5. 8. Results for specific theories 

Before closing this section, we explicitly apply our results to two interesting 3 + 1 
dimensional theories, pure Yang-Mills theory and the A/" = 4 SYM theory. Using the 
single-particle partition functions given in ( |3.20| ), we may write down explicit expressions 
for the transition temperature using z{xh) = 1, for the behavior of the free energy near 
the transition using ( |5.51j ), and for the high temperature behavior using ( |5.18| ). 

For pure (i = 3 + 1, U{N) (or SU{N)) Yang-Mills theory on an of unit radius, the 
single-particle partition functions are given by 

zb(x) = _ , zf{x) = 0. (5.57) 

[L X) 

Solving z{xh) = 1, we find that the Hagedorn temperature is given by xh — 2 — y^, 
Th = -l/ln(2 - VS) - 0.759326. The free energy below the Hagedorn temperature 
may be obtained from the partition function Z in (|3.8|) using ( ^.57|) and recalling that 
X = exp(— 1/T). Just above the phase transition, the free energy is given by 

-^FAr=o(T-^T+) = -0.9877(T-Th)-3.004(T-Th)^-5.980(T-Th)' + O((T-Th)^), 

(5.58) 

while the high temperature free energy is 

Fj^=o{T » Th) ^ -^iV^TVsa, (5.59) 

45 

where Vss = 2ii'^R^ is the volume of the sphere (restoring a general value for the 
radius). 

For the d = 3 + lAf = 4 SYM theory on an of unit radius, the single-particle 
partition functions are 

, , 6x + 12x2-2x^ , , 16xi ,^ 

zb{x) = -3 , zf{x) = -3, 5.60 

(1 — X)"^ (1 — x)"^ 

so from zb{xh) + zf{xh) = 1 we find that the Hagedorn temperature is given byil 
Xh = 7 - 4^3, Th = -l/ln(7 - 4^/3) ~ 0.379663. The free energy below the Hagedorn 
temperature may be obtained from the partition function Z in ( |3.8| ) using (|5.60|) , with 



"^^ Curiously, xh for the TV = 4 theory is exactly the square of xh for the pure Yang-Mills theory, 
and Th in TV = 4 SYM is precisely half of Th in pure Yang-Mills. This arises from the even more 
peculiar fact that the single-particle partition functions are related by zj\f=o{x) = z^=4{x^). 
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X = exp(— 1/T). The free energy immediately above the phase transition temperature is 
given by 

^Fa^=4(T^T+) = -0.9877(T-Th)-4.248(T-Th)^-11.696(T-Th)'+O((T-Th)^), 

(5.61) 

while the high temperature free energy is 

Fj^=4{T > Th) ^ -yiV^TVsa. (5.62) 

6. Phase structure at weak coupling 

In this section we turn to an analysis of the phase transition at non-zero 't Hooft 
couplingii. We will see that an arbitrarily small non-zero coupling qualitatively changes 
the behavior of the transition. In §6.1 we deduce the most general form of Seff{U,X) 
allowed by gauge invariance. In §6.2 we analyze the structure of Seff{U, A) in perturbation 
theory. In §6.3 we study the structure of the deconfinement transition (discussed in detail 
for the free theory in the previous section) at weak coupling. In §6.4 we verify the general 
analysis of §6.3 in an exactly solvable toy model. We end in §6.5 with a description of the 
implications of our results for the micro canonical ensemble. 

6.1. General properties of the effective action 

For any value of the coupling constant, gauge invariance imposes tight constraints on 
the form of the effective action Sef / (a) (|4.6|) for the zero mode of the gauge field; Sef / (a) 
should be invariant under all space-time gauge transformations that 

(1) Are single valued (up to an element of the center of U{N)) on A4 x , 

(2) Preserve the gauge-fixing conditions (|4.1|) and (f4.2|). 

We will restrict attention to gauge transformations U{t) that are independent of the 
position on the compact space. Under such a transformation 

Ao ^ V{t)AoV^{t) - idtVV\ (6.1) 

Some of the key ideas underlying this section and the next one arose in discussions with 
R. Gopakumar. We thank him in particular for emphasizing the qualitative difference between 
the two scenarios depicted in figure 4, and for forcing us to understand the order parameter in a 
clearer fashion. 
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so that 

a V{t)aV\t) - idtVVl (6.2) 

V{t) obeys the condition (2) above (and ( |6.2| ) makes sense) only when the right-hand side 
of (|6.2|) is independent of time. 

Constant (time independent) gauge transformations clearly satisfy the requirements 
(1) and (2) above. Consequently, Seff{a) is invariant under a VaV^ . We may use 
this invariance to diagonalize a. Once this has been done, we consider the further gauge 
transformations 

V{t) = e''^, (6.3) 

where D is a, diagonal matrix, whose eigenvalues are all integral multiples of 2n/(3. a 
transforms under the gauge transformation generated by ( |6.3|) as a a + D. This implies 
that Seff{oi) is invariant under separate shifts of any of the eigenvalues of a by multiples 
of l-K I (5. It follows from these invariances that S^ff is really a function only of tr([/'^) (for 
all n) where U = e*^" (see ( [4.12|) ) is the zero mode holonomy around the time circle. 

Finally, consider V\t) = e 0^ where k is an integer if the gauge group is SU{N), 
and arbitrary for U{N). V{t) obeys the single- valuedness condition (1), as belongs 
to the center of the gauge group. Under the gauge transformation generated by V{t), 
a a + Consequently, S^ff should be invariant under U e^^U, for all integers 
k if the gauge group is SU{N), and for any k for U{N). In the limit N ^ oo the two 
cases coincide, and Seff{U) must be invariant under U — e'^^U for arbitrary 9. Putting 
everything together, we conclude that S^f / may depend on U only in combinations of the 
form 

triU''^) ■ ■ ■tr(C/"'=)tr(C/-'"i---'''=) . (6.4) 

6.2. The general form of the effective action in perturbation theory 

In large perturbation theory, Seff is generated by planar vacuum diagrams ob- 
tained by integrating out the massive modes. In this calculation, a is a background field, 
which appears diagrammatically via external line insertions on the index loops of diagrams. 
Planar diagrams at order A'^"^ include {k + 1) index loops, each of which leads to a trace 
in the final result (which may contain arbitrarily many factors of a) . From the discussion 
of the previous subsection, the only terms with {k + 1) traces allowed by gauge invariance 
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take the form ( |6.4|) . It follows that the most general form of the planar contribution to 
Sef f in perturbation theory is 

m — n—m 

+ C.C.) 

(6.5) 

-|- A ^ ^ -^m,n,p('^) ('^m'^n'^p'^ — m— n— p ^~ C.C.) -|- " • •] , 

m,n,p 

where we have defined Un = tr{U^)/N (note that w* = U-n and uq = 1). This agrees, of 
course, with the quadratic form of the one-loop effective action that we found in sections 
3 and 4. Also, the w„ coincide with the variables pn used in §5 for eigenvalue distributions 
which are symmetric about ^ = 0, so the values of in the free field theory may be read 
off from equation ( p.7| ) . 

6.3. Possible phase structures at weak coupling 

Recall from §5 that, at A = 0, ui is massless at the phase transition temperature Th, 
while the other u^s are all massive. Consequently, ui will continue to be the lightest mode 
in the vicinity of the phase transition also at weak coupling. Thus, for analyzing the phase 
transition at weak coupling, it is useful to obtain an effective action for ui by integrating 
out the u^s with n > 1 @. It follows from (|6.5| ) that the leading large terms in the 
resulting effective action are of the form 

SeffM=N^ ^m?(x,A)|wi|2 + ^^A2"-2s„(x,A)|t/i|2"j , (6.6) 

where ml and the B^s are functions of the temperature and power series in A starting 
(generically) from a A'^ term. 

In particular, to O(A^) we have 

S,ff{ui) = {mjluil^ + b\ui\^) , (6.7) 

where b{x) = i?2(a^,0)A^ is a function of temperature which is generically non-zero. It is 
not difficult to compute -62(0) starting from (|6.5D . The only terms in ( |6.5|) that contribute 
are 

S',ff{U) = {ml\ui\^ + mi|zi2p + A/(w_2W? + W2W-1) + X^A\ui\^) . (6.8) 
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As we are interested in the large N limit, it is sufficient to integrate out these modes classically. 
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At this order, U2 may be integrated out from ( |6.8|) by setting it to be equal to its classical 
value U2 = —IXul/ml il, yielding ( |6?7| ) with b = {A — P/rnDX^. Note that because the 
eigenvalue density (j^^) has to be non-negative everywhere, our expressions for the effective 
action and (|6.7|) are valid only for ui such that |wi| < ^ + 



As in §5, the one- loop contribution to the free energy coming from integrating over 
ui diverges at the temperature Th{X) at which m^(x,A) goes to zero and ui becomes 
massless. For any value of the coupling, at leading order in the distance from this tem- 
perature, we have ml ~ K{Th{X) — T) for some positive constant K. This divergence 
signals a Hagedorn behavior of the single-particle spectrum of the theory with Hagedorn 
temperature Th (computable in perturbation theory), so this behavior persists (at least 
in the microcanonical ensemble) even at non-zero 't Hooft coupling. 

As the saddle point at wi = is unstable for T > Th, the theory described by ( |677| ) 
clearly undergoes a phase transition (to another saddle point) at some T < Th- Whether 
this phase transition occurs at T < Th or T = Th depends on whether the value of b 
(defined in (|6.7|)) at the Hagedorn temperature Th is positive or negative, as we now argue 
in detail 0. The formulas in the remainder of this section are all correct only to leading 
order in b (or in A). 

First, consider the case 6 > 0. For T < Th, ui = 0, corresponding to the uniform 
eigenvalue distribution, is clearly a global minimum of the effective action (|6.7| ). For 
T > Th, however, wi = is unstable and Seff is minimized at = |mip/26. The 

value of the effective action at this minimum is Seff = —A^^|mi 1^/46, which is of order 
(T — ThY , and so the phase transition at T = Th is of second order. As the temperature 
is raised above Th, the eigenvalue distribution smoothly becomes non-uniform until we 
reach \ui\^ = 1/4; this occurs at mf = —6/2. At this point, p{9) vanishes at some 9 and 
we have reached the boundary of the space of eigenvalue distributions and the edge of the 
validity of (|6.7|). As the temperature is further raised the eigenvalue distribution develops 
a gap on the circle, and the theory undergoes a further phase transition similar to the 

Thus setting the eigenvalue distribution to 




(6.9) 



for some a. 



lee [26 1 for a similar discussion in a related context. 
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Figure 2: Plots of Seff{ui) as a function of ui for small positive b, in units 
of N'^b, at several values of mf (from top to bottom) : = (Hagedorn 
temperature and first phase transition temperature), = — |, mf = — | 
(second phase transition temperature), and ml = 

Gross- Witten transition 0. This second phase transition (at T = Th + 27c) is of third 
order. The behavior above this second transition temperature is no longer captured by the 
effective action (|6.6|), and the full action is required to analyze it. 

In the case where b is negative, Seff develops a new local minimum at |wi| = | (the 
boundary of our order parameter space) when mf < \b\/2; note this happens below the 
Hagedorn temperature. When > \b\/4, the free energy at this new local minimum is 
positive, and the saddle point = ^ is disfavored compared to ui = 0. However, when 
we raise the temperature to ml < \b\/4, the free energy at |wi| = | becomes negative, and 
so this saddle dominates over the ui = saddle point. Consequently, at T = Th — 
|wi| jumps discontinuously from zero to ^ and the theory undergoes a first order phase 
transition. When the temperature is further raised, the eigenvalue distribution develops 
a break on the circle and the theory is no longer described by ( [6.61 ). This behavior is 
qualitatively similar to that of the free theory which we analyzed in §5, except that the 
phase transition now happens below the Hagedorn temperature Th- 

If b vanishes exactly then the higher order terms in ( |6.6| ) are required for analyzing 
the phase transition, but generically b will not vanish at the Hagedorn temperature. In 
our discussions above we assumed that b was constant in a range of temperatures near the 
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Figure 3: Plots of Seff{ui) as a function of ui for small negative 6, in units 
of A^^|6|, at several values of ml (from top to bottom) : = m\ = -^j- 
(new phase nucleated), m\ = ■'j- (first order phase transition temperature), 
m\ = -^l^-, and m\ = — 

Hagedorn temperature; this is consistent at weak coupling since the range of temperatures 
which is relevant for the discussions above is of order A^, so the changes in h within this 
range are of higher order in A. 

6.4- Toy model 

The varieties of behavior described in the previous subsection may be illustrated using 
a simple toy example, given by the matrix model 

Z{(3) = j [dU] exp [- (|tr(C7)|2(m2 - 1) + 6|tr(C7)|ViV')] ■ (6.10) 

Based on the previous discussion, we expect this toy model to give a good picture of the 
behavior of the full theory for small A near the phase transition, since the u^s for |n| > 1 
are all small there. Readers who believe our general discussion above may skip directly to 
the next subsection. 

Our toy model may be solved exactly using the methods of §5.4. Since the action 
in ( |6.10| ) depends only on tr(t/) and not on traces of higher powers of U , the possible 
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forms of the saddle point eigenvalue distribution for the integral ( |6.1C1| ) (up to an overall 
shift in 9) are the same as those in (and in our truncated model in §5). The possible 
distributions may be parametrized by a positive number g, and they are given by one of 
the two expressions 



p(^) = ;lfl + -cos(^)) g>2, 
In \ g 



(6.11) 



where for the matrix model ( p.lO|) we have 



g-^ = ui{l-ml~ 2bul) (6.12) 

with 

u^ = ^-^= f dep{e)cos{e). (6.13) 



Note that ui is real since we assumed in ( |6.11|) that the eigenvalue distribution is symmetric 



about 6' = 0. 

We search for simultaneous solutions to (|6.11|) , (|6.12|) and ( |6.13|) . It is easy to check 
that the constant distribution ui = Q {g = oo, p{6) = l/27r) is always a solution; we call 
this solution Phase I. The free energy in Phase I is zero. From the discussion of §5 we 
know that this solution is a local minimum when > 0. 

An additional solution (Phase II) is obtained by using the first distribution in ( |6.11J ) 
(with g > 2), upon solving (|6.12| ) and ( |6.13|) for this case. This gives 

Ml - (1 - ml)ui - 2bul, (6.14) 

with the solution 

«f = -^. (6.15) 

Obviously, this only makes sense if m\/h < 0; consistency of the solution also requires 



2 

mf 



2b 



< \. (6.16) 



so this phase exists and is a local minimum for 6 > and —b/2< ml < (such a phase 
exists also for ml > and 6 < 0, but in that case its free energy is positive and so this 
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phase is sub-dominant compared to Phase I). The free energy F = — Tln(Z) evaluated on 
this saddle point is 

F{T) = (6.17) 

which is quadratic in (T — Th) near the Hagedorn temperature (where oc Th — T). 

To find solutions using the second distribution in ( |6.11| ) (with g < 2), we must once 
again solve ( |6.12|) and ( |6.13|) with the appropriate distribution, and we find 

Given the solution for wi, the free energy evaluated on this saddle point is 
F{T) = N^T 



{ml -l)ui + but + \-^-ln{2il-m)) 



(6.19) 



(the first two terms come from the potential in ( |6.10| ) evaluated on the saddle, while the 
remaining two terms come from the measure). 

The exact solution to the quartic equation (|6.18|) is easily obtained, though it is 



not particularly illuminating. It may be checked that there are no physical solutions for 
rnf > — |, while for mf = — | we find ui = ^. The discussion of the previous subsection 
suggests (and we will demonstrate this below) that this solution represents the phase 
transition point from Phase II to Phase III when 6 > 0, and a phase nucleation point (the 
creation of a new phase at the boundary) when 6 < 0. 

To understand the behavior near this value of , we analyze the solutions to ( |6.18| ) 
in the neighborhood of this special point. For this purpose we introduce a new variable 
y = —ml — I which is linear in 5T, the increment over the phase transition (or nucleation) 
temperature. We also introduce a variable 5 = ^{4:u1 — 1) that measures the deviation of 
Ui from its value at the phase transition (or nucleation) point. In terms of these variables, 
( |6.18|) may be rewritten as 

l+y-d) = l. (6.20) 

The explicit solution for y in terms of S follows immediately from this; this solution may 
be expanded in a power series in d and inverted to obtain 

s = y-C + ^-^y' -(^ + ^ + 4j)y' + oiy% (6.21) 




4^ The following three formulas are only correct for b > 0; for 6 < there is another solution to 
(|6.20| ) which gives the leading contribution p^ . We thank L. Alvarez-Gaume, C. Gomez, H. Liu 
and S. Wadia for pointing this out to us. 
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Using our definition 5 = |(4w,i — 1) we now find 



ui = - \ 1 + -y -^y^ H -T y + C(y 

2 V b 2b^ 265 ^ ' 



(6.22) 



and inserting this into ( |6.19|) we obtain 
F{T) = N'^T 
= N^T 
= N^T 



16 - 4 - 46 + 66^^ 



1 



, 1 „,3 



(6.23) 



In summary, Phase I exists for all values of mi and 6. Its free energy is zero. Phase 
II exists when < —m\/b < 1/2. Its free energy is F = -N^Tmf/4b. Finally, Phase III 
exists whenever there is a solution with ui > 1/2 to ( |6.18| ) and has a free energy given by 
( |6.19| ) (for 6 > this phase exists when —m\/b > 1/2, and for m\ close to —6/2 it has a 
free energy given by ( |6.23| )). Recall that = K{Th — T) where K is positive. 

It is now clear that the dynamics of our toy model near the Hagedorn temperature 
uiKTh) = depends sensitively on the sign of 6 at this temperature. First consider 6 > 0. 
For T < Th {Tn\ > 0), Phase I with wi = is the only saddle point, and it has vanishing free 
energy. As T rises above Th (for small negative ml) a new phase (Phase II) comes into 
existence. Simultaneously Phase I becomes unstable, and the theory executes a second 
order phase transition into Phase II at T = Th- Finally, at a still higher temperature 
(mf = —6/2) Phase II evolves continuously into Phase III; this transition is of third order, 
due to (|6.17|) , ( |6.23|) , and the fact that y in (|6.23|) is linear in ST. As the temperature is 
raised further, the theory remains in Phase III. 

On the other hand, if 6 < 0, Phase I is the only saddle point for ml > —6/2. At 
this temperature (below the Hagedorn temperature) a new phase comes into existence. 
As the temperature is further raised, this new phase splits into two different saddle points 
(Phase II and Phase III). Phase II has positive free energy whenever it exists, and is always 
unstable. Phase III is locally stable; however its free energy is positive in the neighborhood 
of ml = —6/2, and so (at this point) this saddle is sub-dominant compared to Phase I. 
Upon further raising the temperature, F{T) of Phase III becomes negative at some positive 
value of mf, which may be obtained by solving (|6.18| ) and ( |6.19| ) simultaneously for ui and 
mf] note that this happens below the Hagedorn temperature. Finally, as the temperature is 
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raised above the Hagedorn temperature, Phase II disappears and Phase I becomes unstable, 
so the theory remains in Phase III at all higher temperatures. 

In conclusion, the exact solution of the toy model ( |6.10| ) verifies, in complete detail, 
the general results of the previous subsection. 



6. 5. Density of states as a function of energy 

In the previous subsections we discussed the canonical partition function as a function 
of temperature for Yang-Mills theories at weak coupling. We found three qualitatively 
different classes of Z{(3) depending on whether b (see (|6.7|)) is negative, zero or positive. 
In this subsection we will qualitatively describe the corresponding Yang-Mills theories in 
the microcanonical ensemble, seeing how the density of states depends on the energy in 
each of these three different classes of models. 

In order to make contact with the canonical ensemble discussed in previous sub- 
sections, we will find it convenient to characterize the microcanonical density of states 
S{E) = ln(p(i?)) in a rather peculiar way. We will find it convenient to plot the logarithm 
of the effective temperature \n{T(E)) = —\ii{dS{E)/dE) as a function of ln{E) such 
plots were used for very similar purposes in p3[| . 

Log(T) 



Log(T) 



b<0 




Log(E/N-) 



Log(T) 



Log(E/N-=) 




Log(E/N-) 



Figure 4: Plot of log(T) as a function of log(£') (for energies of order A^^, 
in the microcanonical ensemble), for A = and for A > with 6 > and 
b<0. 
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We thank R. Gopakumar for suggesting this to us. 
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The behavior of the density of states (or ^) as a function of temperature may be 
deduced directly from the behavior of F{T) using the relation E = dp{(3F). However, this 
information is incomplete where there are first order phase transitions, since these corre- 
spond to temperatures at which the energy jumps discontinuously in going from one phase 
to another (non-zero latent heat). For the case A = 0, we have seen that the first order 
transition is a Hagedorn transition, so we may conclude that for micro canonical energies 
intermediate between the two phases, the density of states is Hagedorn, corresponding to 
a horizontal line in the graph of ln(T) as a function of ln{E). Thus, the F(T) behavior 
depicted in figure 1 leads to the microcanonical behavior shown on the left in figure 4. For 
small A > 0, the density of states should be a small perturbation of this behavior, and 
we may deduce from the results of §6.3 that the qualitatively different behavior obtained 
for 6 > and 6 < corresponds to whether the flat Hagedorn region in the A = plot 
tilts upwards (6 > 0), resulting in two continuous phase transitions at Th and some higher 
temperature Ti, or downwards (6 < 0), preserving the first order transition (dashed lines 
in figure 4 indicate points not accessible in the canonical ensemble). 

The plots in figure 4 correspond to energies of order A^^, and they are drawn in the 
large limit where the canonical ensemble exhibits sharp transitions and the boundaries 
between phases in the microcanonical ensemble are distinct. To understand in more detail 
the phases available in the microcanonical ensemble at general energies (including energies 
of order one), it is convenient to stretch out the low energy regions of the three diagrams 
for large but finite N, as we do in figures 5, 6, and 7 below .El 

In the case 6 > this plot is displayed in figure 5. As is apparent from the figure, at 
any given temperature the theory has exactly one phase available to it; this is consistent 
with the behavior described by the series of Landau- Ginzburg diagrams in figure 2 and 
with the analysis of the toy model of the previous section. For T < Th the theory is in 
Phase I. At Th, the theory undergoes a second order Hagedorn phase transition, emerging 
into Phase H (the energy after the transition is of order A^^). At a still higher temperature, 
the theory undergoes yet another continuous phase transition into Phase HI. 

Note that Phase I corresponds to ui = 0, Phase H corresponds to the minimum of (for 
instance) the second plot in figure 2, and Phase HI is at the boundary of the configuration 
space of eigenvalue distributions. 

"^^ Strictly speaking, the boundaries between phases are smoothed out when is finite, but for 
large enough these should look sharp as shown. 



55 



-Log(dS/dE) 
= Log(T) 



Phase 111 



Hagedorn 



Phase II 



ifhase I 



Log(E) 

Figure 5: Plot of log(T) as a function of log(i?) (in the microcanonical 
ensemble) when 6 > 0. 



-Log(dS/dE) 
= Log(T) 



Hagedorn 




Log(E) 

Figure 6: Plot of log(T) as a function of log(i?) (in the microcanonical 
ensemble) in the free Yang-Mills theory. 

The density of states of the free theory may be obtained from the limit 6 ^ of the 
previous discussion; the relevant plot takes the form shown in figure 6. As is apparent from 
this figure, the free theory makes a phase transition (of first order) directly from Phase I 
to Phase III at the Hagedorn temperature, in agreement with the analysis of §5. 

When 6 < 0, we believei2l that the density of states is characterized by the plot 



As noted above, when a theory undergoes a first order transition, Z{P) is insufficient to 
reproduce the full density of states of the theory. In order to obtain figure 7 we use one additional 



56 



-Log(S'(E)) 
= Log(T) 

Log(T_2) 
Log(T_l) 



1 

p/sel 


Hagedom 




Hiase III 






Phasejl / 
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Log(E) Log(E_l) Log(E_2) 

Figure 7: Plot of log(T) as a function of log(£^) (in the microcanonical 
ensemble) when b < 0. 

displayed in figure 7. The thermodynamics induced by the density of states plotted in 
figure 7 is consistent with the behavior described by the Landau Ginzburg plots in figure 
3. In particular, the five graphs shown in figure 3 (from top to bottom) should be taken to 
represent SeffM for T < Ti, T = T^, T = T2, T2 < T < Th and T > Th, respectively. 
Phase I lies at tti = in figure 3. Phase II is the unstable maximum in the 3rd and 4th 
graphs in figure 3, and Phase III is the saddle point at the boundary ui = ^. Note that, 
according to both figures 7 and 3, Phase III exists only for T > Ti, and Phase II exists 
only for Ti < T < Th- 

In ending this section we would like to emphasize that adding interactions has two 
qualitatively different effects on the spectrum of weakly coupled Yang-Mills theories. First, 
single-trace states involving a number of fields which is much smaller than N pick up 
corrections to their energy, starting at first order in A (for conformal field theories on 
S'^~^ these corrections are equivalent to the anomalous dimensions of the single-trace 
operators). These corrections result in a renormalization of mf(T), and hence of the 
Hagedorn temperature. However, the physics of the phase transition (and, in particular, 
the crucial sign of the coefficient b) is governed by a different effect, which is the deviation 

assumption, that the unstable saddle points appearing in the Landau Ginzburg diagram of figure 
3 (see also the toy model of the previous subsection) may be interpreted as an unstable phase (of 
negative specific heat) at the corresponding temperature, in the microcanonical ensemble. 
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from Hagedorn-type behavior in states which involve a number of fields of order A^^. 
This effect appears even in the free Yang-Mills theory, and it appears to be modified in a 
qualitatively important fashion by coupling constant effects at (9(A^), as shown in figure 
4. These modifications determine the sign of b and the nature of the phase transition at 
finite coupling. 



7. Extrapolation to strong coupling and the dual description 

In this section we will present possible extrapolations of the results of section 6 to 
strong coupling. We will also discuss the possible interpretation of these results in a 
stringy dual description. Much of this section is rather speculative; we will present some 
questions and outline possible answers. 

7.1. Possible phase diagrams for large N Yang-Mills theories 

Consider the deconfining transition of a d- dimensional confining large Yang-Mills 
theory on a compact space of size R, at large values of RAqcd oo. As we reviewed in 
§2.2, this transition may be either of first order (as seems to be the case for pure Yang-Mills 
theory) or of second order (in which case, for d — 4 it must be in the universality class of 
thed = 3 XY model). 

On the other hand, as we saw in the previous section, at weak coupling (small RAqcd) 
the transition may be of either first or second order, depending on the (perturbatively 
computable) sign of b for the theory in question. For each of the four possibilities for the 
order of the phase transitions at large and small RAqcd we have sketched the simplest 
possible interpolation between the two limits, and the corresponding phase diagrams are 
displayed in figures 8-11. Recall that Phase I is a confined phase, characterized by a uniform 
distribution of eigenvalues for the holonomy matrix On the other hand Phases II and 



When the number of fields is of order ^/N or larger, single-trace and multi-trace states mix, 
and one cannot distinguish between them. This was highlighted recently in the study of the effects 



of non-planar graphs in the Yang-Mills/plane wave duality [44|. 

In confining theories away from weak coupling, when the size of the compact space becomes 
large compared to the scale set by Aqcd, it is not as useful to focus on the zero mode of the 
gauge field, since other modes become light compared to Aqcd- However, while Seff{U) may be 
less relevant in these more general cases, we can consider instead the effective action Seff{U{y)), 
where U (y) is a spatially-dependent unitary matrix given by the Wilson loop around the thermal 



58 



Hagedorn Temperature 




PHASE I 



RAoCD 



Figure 8: The simplest phase diagram for a compactified confining theory 
with negative b and a first order phase transition at R — oo . Only solid lines 
represent phase transitions. 

A 




Hagedorn 

PHASE I 



R AocD 



Figure 9: The simplest phase diagram for a compactified confining theory 
with positive b and a first order phase transition at i? = oo . Only solid lines 
represent phase transitions. 



circle at a point y (see [45| for a related discussion). For large theories, we expect that the 
saddle point configurations for every temperature (in all phases) should be spatially homogeneous, 
so that the corresponding U (y) is constant. Thus, we may still characterize the various phases by 
a single eigenvalue distribution corresponding to this constant U. 
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PHASE I 



R AocD 

Figure 10: The simplest phase diagram for a compactified confining theory 
with negative b and a second order phase transition at = oo . Only solid 
lines represent phase transitions. 




Hagedora 

PHASE I 



R AocD 

Figure 1 1 : The simplest phase diagram for a compactified confining theory 
with positive b and a second order phase transition at i? = oo . 



Ill are both deconfined phases, which are distinguished by the fact that in Phase II the 
eigenvalue distribution is non-vanishing on the whole unit circle, while it vanishes outside 

60 



an interval in Phase III. 

Note that in each of the figures 8-11 all phase transition temperatures at strong 
coupling are of 0{Kqcd)] consequently, in this regime the dimensionless phase transition 
temperature TR increases linearly with the coupling parameter RAqcd- 

Since lattice data indicates that 3 + 1-dimensional pure Yang-Mills theory undergoes 
a single first order deconfining phase transition at large A^, the simplest possibilities for 
the phase diagram of the 3 -|- 1-dimensional large pure Yang-Mills theory on a compact 
space such as are depicted in either figure 8 or figure 9; as above, the value of b for this 
theory (which we will report on in [|10[) will distinguish between these two options. 

Above we drew the simplest possible interpolations between the weak coupling and 
strong coupling behaviors. These simple interpolations are consistent with the values of 
all of the order parameters we discussed at weak and strong coupling, but several other 
rather natural phase diagrams may also be drawn. For instance, figures 8-11 do not apply 
to any theory that undergoes a phase transition as a function of the coupling at zero 
temperature. Such a phase transition is certainly possible, and could be inserted into any 
of the diagrams of figures 8-11 as a vertical line that divides Phase I into two regions, 
distinguished, presumably, by an order parameter unrelated to the Polyakov loop. A 
case where we know that this happens is when the infinite volume theory exhibits global 
symmetry breaking due to quantum effects (for example, the d = 4 = 1 SU{N) SYM 
theory, in which a chiral '^2N symmetry is spontaneously broken to ^2 at infinite volume). 
Since we see no sign of such a symmetry breaking in our analysis at weak couplingEl, it 
seems clear that in such cases there are additional phase boundaries separating weak 
coupling and strong coupling regimes at low temperature. In general, we do not expect 
any such boundaries for asymptotically free gauge theories at high temperatures, since in 
this limit, the theories may be studied perturbatively for all values of RAqcd- 

We now turn to the large SU{N) A/" = 4 SYM theory on an of unit radius and at 
strong 't Hooft coupling A. Using the AdS/CFT correspondence, this theory is equivalent 
to type IIB string theory at weak string coupling on an AdS^ x 5"^ space (using global 



Alternatively, as mentioned for instance in |46], they can be distinguished by the behavior 
of the Polyakov loops for higher representations of the gauge group, such as the fc'th product of 
fundamental representations for large k. 

Of course, as in our discussion above of the breaking of the ^at symmetry associated with 
confinement, we do not expect to see symmetry breaking at finite A'^ and finite volume, but we 
expect to see a sum over different configurations in which the symmetry is broken. 
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coordinates for AdS) with a large radius of curvature. The thermodynamical analysis 
of this theory shows that it undergoes a single first order phase transition as a function 
of temperature, which occurs at Thp = ^ — 0.477465 In addition, since the 

string coupling is small, we have a large range of energies with a Hagedorn behavior of 
the spectrum, with the Hagedorn temperature scaling as Th oc oc A^/^ and going 

to infinity in the limit of strong 't Hooft coupling. 

On the other hand, as we have described in detail above, at A = this theory also 
undergoes a single phase transition as a function of temperature. This phase transition, at 
the temperature Th = — l/ln(7 — 4v^) — 0.379663, coincides with the Hagedorn transition 
and it is weakly of first order. And, as discussed in §6, the behavior of the theory at small 
A depends crucially on the sign of b (defined in ( |6.7|) ). If b turns out to be negative, the 
theory undergoes a single first order phase transition below the Hagedorn temperature. If 
6 > 0, the theory undergoes two continuous phase transitions, the first of which is at the 
Hagedorn temperature. 



TR 



PHASE III 



Hagedorn Tehiperalure 



PHASE I 



Figure 12: Conjectured phase diagram for the A/" = 4 SYM theory on a 
sphere if b is negative. Only solid lines represent phase transitions. 



As we have described above, the value of b for this theory is determined by a set of 



two-loop and three-loop graphs. We will present the result of this computation in ||10|| . 
If b turns out to be negative, the behaviors at weak and strong coupling are similar and 
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Hagedorn , 

PHASE in 




Hagedorn 



PHASE I 



Figure 13: Conjectured phase diagram for the A/" = 4 SYM theory on a 
sphere if h is positive. Only sohd hnes represent phase transitions. 

the phase diagram of the theory has a natural interpolation (shown in figure 12) for all A 
(though more complicated phase diagrams are also possible). 

On the other hand, if h turns out to be positive, the simplest possible phase diagram 
takes the form shown in figure 13. Note the existence of a tri-critical point at a special 
value of the 't Hooft coupling, at which the deconfinement phase transition changes from 
being second order to being first order, and the existence of a new phase at weak coupling 
and intermediate temperatures, il 

Figures 12 and 13 are very similar to figures 8 and 9, with one important difference : 
the phase transition temperature, measured in units of 1/-R, has a finite strong coupling 
limit for the A/" = 4 theory, but increases without bound for confining theories. 

Note that the interaction potential between a quark and an anti-quark, at a distance 
L <^ i?, is expected to behave rather differently in the various phases described above. In 
the case of confining gauge theories at low temperatures, this potential will be Coulomb- 
like at weak coupling (i?A ^ 1) but linear (string-like) for KK ^ 1, as long as L ^ 1/A. 



Figures 12 and 13 have been drawn under the assumption that the large N J\f = A theory 
in flat space does not undergo a phase transition as a function of the coupling. Such a transition 
could perhaps occur, and then it would be depicted by a vertical line that would divide Phase 
III into two diff'erent regions in flgures 12 and 13. We thank R. Gopakumar for emphasizing this 
possibility to us. Similarly, it is possible in principle that Phase I would be divided into difi^erent 
regions. 
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However, in Phase III (and presumably in Phase II as well) the quark- anti- quark force is 
exponentially screened by the intervening plasma. In the case of the A/" = 4 SYM theory 
the quark- anti- quark force is Coulomb-like at every coupling in the 'confining' Phase I, 
but decays exponentially at sufficiently large distances (compared to a scale set by the 
temperature) in the deconfined phases. 

1.2. Dual interpretation of the Af = 4 SYM thermodynamics at strong coupling 

We noted in §2 that it should be possible to understand the possible phase diagrams 
presented earlier in this section in terms of a dual string theory description of the relevant 
gauge theories. Unfortunately, among the 3-1-1 dimensional gauge theories with adjoint 
fields, the only one whose string dual is known is the A/" = 4 supersymmetric Yang-Mills 
theory (and theories related to it by renormalization group fiows) , and this dual is mostly 
understood only in the limit of strong 't Hooft coupling. In this subsection we will argue 
that, at least in this one case, the thermodynamics of the string theory dual fits rather 
nicely with the picture presented in this paper. Previous discussions of the thermodynamics 
of 7\A = 4 SYM as a function of the coupling, which are consistent with ours, appear in 
14^ , and a detailed discussion of the transitions discussed in this subsection appears in 

US®- 



The string dual to A/" = 4 SYM on 5"^ x IR (where the is taken to be of unit radius) 
is type IIB string theory on AdS^ x S^, whose metric may be written as 

ds^ = Rl{- cosh^ p dr^ + dp^ + sinh^ p dflj + dflf) . (7. 1) 

Under this duality the Hamiltonian of the Yang-Mills theory is identified with the generator 
of global time translations dr in the geometry ( |7.1| ). The energy E of the gauge theory 
on a sphere of unit radius is related to the proper energy in string theory at p = by 
Eprop = E/Rq. The radius Rq of the AdS^ x space is related to the 't Hooft coupling 
A and the inverse string tension a' by Rq ~ X^Vo/ (we will ignore all numbers of order 
unity in the qualitative discussion of this subsection). 

For large Ro/\/c7 (large 't Hooft coupling) the density of states p{E) of type IIB 
string theory on ( |7.1| ) has four distinct regimes (see and references therein). The only 



states in the spectrum with proper energy!^ smaller than the string scale, E <^ A^/*, are 



Recall that states of finite energy on ( [7.1D are all localized about p = 0. AdS^ space behaves 
effectively like a four dimensional box of physical radius Rq, as is apparent from the fact that the 
spectrum of dr is discrete with discretization step unity. 
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ten dimensional supergravitons0 whose entropy scales as S{E) oc E^^^^. For E ^ X^^^, 
excited string states are added to the spectrum; the contribution of these states to the 
entropy is 5" ~ EpropVa' ~ X~^/'^E (see §2.1), and this contribution dominates over the 
graviton gas for E ^ A^/^. When the proper energy exceeds the Planck mass (namely, 
E ^ nisRo/ gV^ = N^/"^) small ten dimensional Schwarzschild black holes are also added 
to the spectrum (note that these energies are inaccessible in the strict N ^ oo limit). 
The entropy of such black holes is proportional to ~ {IpE / R^f'^'^ = (E'/A^^/^)^^^; they 
have negative specific heat and positive free energy. This entropy of black holes dominates 
over that of the Hagedorn strings for E ^ A^^/A^/^. Finally, at an energy Ei ~ A^^ 
(an energy at which the radius of these Schwarzschild black holes becomes comparable 
to the AdS radius Rq) the black hole horizon covers the whole 5"^, and the specific heat 
of these black holes becomes positive. For E > Ei these black holes are referred to as 
AiiS'-Schwarzschild or big black holes. At a higher energy E2 (also of 0{N'^)) the free 
energy of these black holes becomes negative. Finally, for E ^ A^^, the entropy of these 
black holes is 5" oc N^^'^E^/^, so the thermodynamics of these big black holes resembles 
that of a four dimensional conformal field theory. 

-Log(S'(E)) 
= Log(T) 



Log(T_2) 
Log(T_l) 



Log(E) Log(E_l) Log(E_2) 

Figure 14: log(T) as a function of log(i?) (in the microcanonical ensemble) 
for all the 'phases' of type IIB string theory on AdS^ x S^, when Rq ^ y/c/ . 

In terms of the gauge theory, the only states in this regime are those created by the chiral 
primary operators, their products and their descendants. 
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For comparison with the canonical ensemble it is most convenient (see §6.5) to display 
the density of states described in the previous paragraph as a plot (see figure 14) of 
ln(T(i?)) = — ln{S'{E)) |^^. Quite remarkably, figure 14 is identical in its general form to 
figure 7, even though figure 7 was derived for a general gauge theory (with 6 < 0) at weak 
coupling while figure 14 applies to the A/" = 4 Yang Mills theory at strong coupling, and 
has been derived using its stringy dual. 

To complete this subsection we now use figure 14 to discuss the behavior of the N" = 4 
SYM theory at strong coupling in the canonical ensemble (see §6.5 for a very similar 
discussion). For T < Ti(= ^) the only available saddle point is the thermal gas of 
gravitons. Over the temperature range Ti < T < Th — the theory has three saddle 
points to choose from. One of these saddle points (small black holes) is unstable, so 
it cannot appear in the canonical ensemble. The other two saddle points are stable and 
compete with each other. The gas of gravitons has lower free energy for Ti < T < T2(= ^), 
while the big black hole dominates the canonical ensemble for T2 < T < Th ■ Consequently, 
the theory should undergo a first order phase transition at T = T2, which (as we mentioned 



above) is indeed the case [|13| , pT[| . T = Th is the Hagedorn temperature at which the 
graviton gas saddle point stops existing. For T > Th the generic state of the theory 
continues to be a big black hole, which is now the only available phase. 



7.3. Deconfinement and black holes 

We do not, as yet, understand the string dual of weakly coupled gauge theories, such 
as pure Yang-Mills or A/" = 1 supersymmetric Yang-Mills on a small S^, or even the A/" = 4 
SYM theory at weak 't Hooft coupling. Nonetheless, each of these theories is continuously 
related to the Af = 4 theory at strong coupling, so the AdS/CFT correspondence provides 
a demonstration (in the physicist's sense of the term) of the existence of string duals for 
all these gauge theories. Furthermore, these string duals may be expected to qualitatively 
resemble type IIB string theory on AdS^ x S"^, in the sense that they should all include 
at least one additional "holographic" dimension beyond the dimensions of the Yang-Mills 
theory. Consequently, even though we understand little about the details of the string 
duals of arbitrary gauge theories, we will use intuition from the previous subsection to 
boldly speculate on the stringy dual interpretation of the phase diagrams in figures 8-13. 
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As we reviewed above, the deconfined phase of the strongly coupled A/" = 4 theory 
admits a dual description in terms of black holes il. We would now like to argue that this 
is a general feature (see for related comments): large deconfined phases (Phases II or 
III in figures 8-13) should always be associated with black holes ii. Our simple argument 
follows directly from ideas presented in |]Tl[ . 

In order to measure the Polyakov loop we need to put an external quark on a trajectory 
that wraps around the circle of the time direction. Since in the large N gauge/string 
correspondence quarks appear at the boundaries of the worldsheet, we expect the Polyakov 
loop in string theory to be realized as the partition function (in Euclidean space) for a 
string whose boundary wraps around the time-like loop, and this is indeed the case in 
the AdS/CFT correspondence [56,5^11. This partition function clearly vanishes if the 
Euclidean time circle is non-contractible0, but it is generically non-zero (when suitably 
defined, as in our discussion of §5.7 and the analogous discussions in |jTl|,^) if the time 
circle is contractible. In all the examples that we know of, a static Euclidean time cycle is 
contractible in general relativity only when the corresponding Minkowskian solution has 
a horizon. Outside the range of validity of general relativity it would seem reasonable 
to define black holes by this requirement. Putting all this together, we conclude that 
suitably defined Polyakov loops are non-zero only in a black hole phase. Thus, we generally 
expect that the string theory (at finite temperature) in Phases II and III will be in a 
generalized black hole background, while the theory in Phase I (where the time cycle is 



The connection between deconfinement transitions and black holes has also been verified in 
examples of confining field theories in fiat space, such as the duality between a cascading d = A 
SU{N + M) X SU{N) TV = 1 supersymmetric gauge theory and type TIB string theory on a 
manifold which includes the resolved conifold pT|,|52|. Similar results for the Maldacena-Nunez 



background [ p3| were obtained in [54|, and for other theories in (for example) [ p5[ . 

By this we mean that if the background has a geometrical interpretation it would resemble 
a black hole; we generalize this to other backgrounds by calling such phases "black holes" . 

In this case the boundary of the string is actually not precisely the Polyakov loop, as it 
contains also couplings to scalar fields. We expect that these additional couplings will not change 
the qualitative features of the order parameter, and that these "generalized Polyakov loops" will 
vanish if and only if the usual Polyakov loop vanishes. Computations of such "generalized Polyakov 
loops" in the AdS/CFT correspondence first appeared in [11,|58|,59|]. 

Again, this is true when we have a good geometrical description of the background, and 
otherwise we take the vanishing of the Polyakov loop to be the definition of a non-contractible 
circle. 
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non-contractible) should be in a background similar to the naive gauge theory background 
(in which the time direction is simply compactified). 

A similar conclusion follows by analyzing the other order parameter for confinement, 
lim.N F{T)/N'^ . In string theory this order parameter maps to the sphere partition 
function, which naively always vanishes. However, this is not true even in the low-energy 



general relativity approximation for some types of unbounded spaces [^|TT|, where the 
naive classical action diverges and one needs a subtraction procedure to define it. This can 
lead to a non-zero result for the relative classical partition function between two spaces 
with the same asymptotics but a different interior. This again suggests that deconfined 
phases must involve different space-times than the confined phase (which are the same 
asymptotically but differ in the interior), which have different classical actions. The dis- 
cussion of the previous paragraph suggests that these spaces should be generalized black 
holes. 



7.4- Dual description at a general point in the phase diagram 

At any fixed value of the coupling in each of the diagrams of figures 8-13, the com- 
pactified large gauge theory either 

(1) Undergoes a single first order phase transition as a function of the temperature away 
from a tri-critical point. 

(2) Undergoes a single first order phase transition as a function of the temperature at a 
tri-critical point. 

(3) Undergoes two successive continuous phase transitions as a function of the tempera- 
ture. 

It is possible (see below for a caveat) that the density of states takes the form shown in 
figures 7, 6 and 5, respectively, in every realization of the three cases listed above. We 
have already argued that this statement is true at weak coupling (see §6.5), and for the 
A/" = 4 SYM theory at strong coupling. 

If this picture turns out to be correct, it would be tempting to further speculate that 
the dual description of Phase III always involves a stable 'big black hole'. Following this 
train of thought, one would also be led to speculate (following §7.2) that Phase II has a 
dual description in terms of a different type of strange 'black holes' (see §7.3 for our usage 
of the term black holes). It would then follow that the discovery of theories of type (3) 
(theories with a second order large N deconfinement phase transition) would imply the 
existence of a new class of stable 'strange black holes' in a class of string theories. This 
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would certainly be intriguing; it would be very interesting to understand the string theory 
interpretation of the distinction between the two classes of 'black holes' (Phase II and 
Phase III). 

Before concluding this section, we should note that the different theories and phases 
which we grouped together in our discussion above do seem to have some qualitative 
differences, despite sharing some similar features. Let us first study the behavior of a 
confining gauge theory (like d = 4 pure Yang-Mills theory) as a function of coupling 
in more detail. On a very large compact space (in units set by Aqcd) the pure Yang- 
Mills theory behaves locally just like the infinite volume theory, and we expect that the 
deconfinement transition proceeds (as energy is added to the system at the transition 
temperature) by nucleation of bubbles of the deconfined phase, which grow and eventually 
coalesce to cover the whole space. In particular, in this strong coupling limit there should 
be stable configurations at the transition temperature characterized by the coexistence of 
confined and deconfined phases covering different regions of the compact space. S In such 
mixed configurations we would expect the size of every phase bubble to be at least of order 
Aqqj^; in particular, mixed phases may be expected to be absent on a compact space that 
is small in QCD units, namely at weak coupling AqcdR ^1.1^ This is consistent with 
the analysis in our paper; recall that the order parameter we used for the phase transition 
is (roughly) the constant mode of {jjtr{U{y))) on the compact space, and that at weak 
coupling modes that describe inhomogeneous configurations of tr{U{y)) are massive with a 
mass of the same order as the phase transition temperature, and may safely be integrated 
out, as we have done. 

The situation is completely different in the case of the d = 4 A/" = 4 SYM theory. In 
this theory the existence of a phase transition relies upon the compactness of the space at 
all values of the coupling, and the behavior at weak coupling and at strong coupling seem 
similar. Recall that, in the strongly coupled A/" = 4 theory, the dual description of the high 

Note that in our discussion above of microcanonical phase diagrams (and, indeed, in our 
whole analysis which was based on the effective action Seff{U)), we were implicitly discussing 
only homogeneous phases. It seems that at least in some cases non-homogeneous configurations 
dominate at strong coupling (in the microcanonical ensemble) over the configurations we described 
above. 

The fact that the qualitative behavior of Yang-Mills theories changes in this respect in the 
transition from weak coupling to strong coupling may indicate that the naive interpolation we 
suggested between these two regimes is too optimistic. 
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temperature phase is a single large black hole. Since there would seem to be no such thing 
as half of a black hole, the high temperature phase cannot coexist with a low temperature 
phase at the transition temperature at strong coupling. 

The contrasting behaviors of the systems described in the previous two paragraphs 
suggest that the "large black hole" phase (Phase III) of strongly coupled pure Yang- 
Mills theory may be significantly different from the large AdS black hole. It would be 
interesting to understand what sort of bulk description the coexistence of confined and 
deconfined phases in confining theories such as the pure Yang-Mills theory could have. 
One is reminded of the appearance of non-homogeneous configurations in the context of 



the Gregory-Laflamme transition [61 



8. Discussion and future directions 

In this paper we have analyzed the thermodynamics of weakly coupled large gauge 
theories compactified on a sphere of radius R, or any other compact manifold on which the 
theory has no zero modes. Our analysis applies both to conformal gauge theories at small 
values of a tunable coupling constant (such as the Af = 4 SYM theory) and to confining 
theories (such as pure or A/" = 1 Yang-Mills theory) with RAqcd ^ 1- We have shown that 
in the microcanonical ensemble these theories exhibit an exponential (Hagedorn) density 
of states, which is cut off at an energy of order E ~ A^ (this was previously shown 
for free gauge theories in [§,|5[). The Hagedorn temperature at zero coupling is easily 
determined from the field content of the theory, and the corrections to it can be computed 
in perturbation theory. We have demonstrated that these theories undergo deconfinement 
phase transitions, with different possible phase diagrams discussed in §6. 

Our analysis has several points of interest. First, putting a confining gauge theory 
on a compact space introduces a new dimensionless parameter into the game; varying this 
parameter may continuously deform the flat space deconfinement transition into a regime 
where the transition (which remains non-trivial) may be reliably studied in perturbation 
theoryi^. The analysis of this weakly coupled problem leads to several insights. For 
instance, we have argued that if the large A deconfining phase transition is of second order 
then it must be Hagedorn-like (as first shown from different arguments in IQ), and it must 

This is somewhat analogous to the continuous relation between confining phases and 
perturbatively-accessible Higgs phases which can occur [ |62| in theories with fields in the funda- 
mental representation upon changing coupling constants or vacuum expectation values of fields. 
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be followed by a second phase transition at a higher temperature (where the eigenvalue 
distribution corresponding to a holonomy around the thermal circle develops a gap). This 
suggests that a second order deconfining transition in any large gauge theory implies the 
existence of a previously unsuspected intermediate temperature phase. Our analysis of the 
weak coupling deconfinement transition led to us suggest the four simplest possible phase 
diagrams (figures 8-11) for confining large gauge theories; which diagram is actually 
implemented depends on the details of the theory, including the sign of a perturbatively 
computable number. More speculatively, we have also suggested a stringy interpretation 
(in terms of nucleation of black holes in the dual string theory) of this deconfinement 
transition. 

Turning now to the A/" = 4 supersymmetric Yang-Mills theory, our results imply that 
the thermal phase transition in this theory at weak coupling is either of first or second order, 
depending on the sign of a coefficient that we are now in the process of computing. If the 
transition turns out to be of first order, the simplest possible conclusion would be that the 
Hagedorn spectrum of string theory on AdS^ x does not dominate the thermodynamics 
of A/" = 4 SYM at any temperature and at any non-zero value of the coupling - this was 



called "Hagedorn censorship" in [^. In such a case our computations would enable us 
to study an unstable saddle point representing substringy unstable ('Schwarzschild') black 
holes in this background. On the other hand, if the weak coupling transition turns out to 
be of second order, that would suggest the existence of a tri-critical point at some finite 
't Hooft coupling (probably of order one) in this supersymmetric Yang-Mills theory, as 
well as the existence of a previously unsuspected intermediate temperature phase in this 
theory at weak coupling. This new intermediate temperature phase should then have a 
dual description in terms of a new set of bulk objects: mysterious new stable black holes. 

The picture we present in this paper supports the view that 't Hooft's relation between 
large gauge theories and string theories [jl| may extend also to weakly coupled gauge 
theories. The naive picture of this relation is that at strong coupling the sum over pla- 
nar Feynman diagrams is dominated by graphs with arbitrarily many interactions, which 
become dense and close up the holes in the worldsheet of the Feynman diagram to form 
a continuous 1 -f 1 dimensional field theory on a spherical Riemann surface. This view 
has already been seriously challenged by the AdS / CFT correspondence, which establishes 
a relationship between string theories and gauge theories on a sphere at all values of the 
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gauge couplingL^, and by the consequent derivation of dualities between (topological) open 
and closed string theories |170| . [71|] . Our findings provide further evidence for the existence 
of a stringy dual of weakly coupled gauge theories. We have found that weakly coupled 
gauge theories on compact manifolds share at least one qualitative feature of string theory 
in flat space, namely a string-like spectrum, in a regime in which the Feynman diagrams 
do not seem to look like continuous Riemann surfaces. It would be very interesting to un- 
derstand more directly how and why weakly coupled gauge theories manage to rearrange 
themselves as string theories (see |l7^ - |7^ for some attempts in this direction). 

It is an interesting challenge to identify the string theories which are dual to various 
weakly coupled (or free) gauge theories. For the free gauge theories our analysis of §3 gives 
us an explicit formula ( |3.6| ) for the spectrum of (free) single-string states, and perhaps in 
some cases this may be enough to reconstruct the corresponding worldsheet theory. By 
expanding our result for the partition function in a power series in we can also extract 

information about string coupling corrections in the putative dual string theory (though 
it may turn out that real string interactions always involve the Yang-Mills interaction, as 
in the study of string interactions on plane waves |[44||). 



The results of the computations in this paper (and our upcoming work [|T0[) may 
turn out to have independent interest. For conformal theories such as the (i = 4 A/" = 4 
supersymmetric Yang-Mills theory, the partition function computed in this paper and in 
||10|| encodes, in principle, the scaling dimension of every operator in the theory, to the order 
of computation (in A) il. The information about these anomalous dimensions is packaged 
in an interesting fashion in Se/f and may lead to insights, perhaps in combination with 
recent speculations on integrability and on the existence of a large Yangian symmetry in 
the AT = 4 SYM theory fT^. 

The techniques of this paper may flnd application to several theories not explicitly 
considered in this paper. It would be interesting to compare our results to known results 
in strongly coupled two dimensional gauge theories pO| , pT|p2| . On another note, recall that 
large N gauge theories in d = 1 with a single scalar field are believed to be holographically 
dual to the c = 1 and c" = 1 string theories. Our analysis of free one-matrix positive-sign 



Although the worldsheet theories dual to perturbative gauge theories are strongly coupled, 
as a consequence of large background curvatures, they presumably exist, as they may be ob- 
tained from the well-understood duals to strongly coupled gauge theories by tuning worldsheet 
parameters. See p^ - p^ ] for recent attempts to understand this weakly coupled limit. 

This is true at least in the absence of first order phase transitions. 



72 



harmonic oscillators at zero coupling shows that they have neither phase transitions nor a 
Hagedorn-like spectrum; however, the matrix models dual to the c = 1 and c = 1 strings 
are 'negative mass' harmonic oscillators, and it is possible that they somehow undergo a 
'deconfining transition' of the sort studied in this paper, described by some effective action 
Sf,ff{ui), and that their high energy behavior is dominated by two dimensional black holes 
(whose spectrum is Hagedorn-like). 1^1 Another interesting example is provided by the 't 
Hooft limit of the maximally supersymmetric plane-wave deformation of Matrix theory, 
argued in to be dual to a little string theory compactified on S^. For a small sphere, 
this theory is weakly coupled, and one may study the little string theory thermodynamics 
explicitly in this limit using methods similar to this paper. This analysis is the subject of a 
paper to appear shortly. Another interesting relation between the thermodynamics of 
string theory, black holes and the thermodynamics of compactified gauge theories appears 
in the context of Matrix theory (see, for instance, |7^). It would be interesting to try to 
relate our discussion to the results for toroidal compactifications which are relevant there. 
When we compactify a gauge theory on a space with a non-zero fundamental group, we 
necessarily have additional zero modes coming from the non-trivial Wilson loops. It would 
be interesting to generalize our results to include such additional degrees of freedom. Some 
recent results on the deconfinement transition of the large gauge theory on a torus in 



the strong coupling limit appear in ||79|| . 

We should also note that the connection between Hagedorn transitions, black hole 
formation and some sort of 'deconfinement' transition in a dual field theory seems more 
general than the specific example of gauge theories. For instance, string theory on AdSs is 
known to undergo a phase transition as a function of temperature both at 'weak coupling' 
Il80[] (at the orbifold point) and at strong coupling (in a geometric phase for the bulk). At 
weak coupling this transition (between a single "long string" phase and a multiple "short 
string" phase) is Hagedorn-like ||81|], while at strong coupling it is a first order Hawking 
Page transition of the sort reviewed in this paper. It would be interesting to analyze 
the extrapolation between these two behaviors. 

Finally, several technical and conceptual issues remain to be addressed even within 
the direct line of attack of our paper. The coefficient b should certainly be computed 



for several gauge theories [10|. The analysis of this paper can be extended to search for 



Several people, including A. Adams, P. Ho, G. Mandal, and S. Wadia, have suggested a 
variant of this idea. 
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interesting features in generalized partition functions t3. And, last but not least, the dual 
interpretation of deconfining phase transitions and their intriguing connection with black 
holes certainly deserves further study. 
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Appendix A. Properties of group characters 

In applications of group theory, it is often necessary to determine the set of irreducible 
representations obtained in the tensor product of some collection of other representations 
Ri. Characters provide a powerful tool to achieve this. Given a group G, we may associate 
with any representation Ri a character '■ G — defined such that XrAU) is equal to 
the trace of the group element U in the representation Ri. From this definition, it follows 
that the characters for sums and products of representations are given by 

XRi®R2 = XRi +XJ?2 5 XRi^R^ = XRi ■ XR2- (A-1) 

The utility of the characters in decomposing tensor products follows from the orthogonality 
of characters for irreducible representations Rj : 

j[dU]x*^,^{U)xRiSU) = ^R{Ri^ (A.2) 

where [dU] is the invariant (Haar) measure on the group manifold normalized so that 
J[dU] = 1. Thus, the number of irreducible representations of type R^ in the tensor 
product of the representations . . . , i?„ is given by 

/n 
[dU]xRi{U)llxRAU). (A.3) 
i=i 

For the special case of the trivial (singlet) representation, we have Xsinglet{U) = 1, so that 
the number of singlets in the product of representations Ri is simply 

nsinglet= j [dU]Y[xR,{U)- (A.4) 

Finally, for the applications of this paper, we require the character formulae for the sym- 
metrized and antisymmetrized products of n identical representations R. To obtain these. 
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let Uji be the matrix representation of the group element U in the representation R. Then 
the trace of the matrix representation of U in the (anti)symmetrized tensor product of n 
copies of R is 

iUR)i^.---iUn)Z]^, (A.5) 

where the [■ ■ ■]± indicates symmetrization or anti-symmetrization of the indices with unit 
weight. This expression is exactly the term in the expansion of the integral 

Gmt) = / [#±]e-^-^±±^^±^«^± (A.6) 

in powers of t, where (/)_(- and cp- are complex bosonic or fermionic variables, respectively, 
in the representation R. Thus, the result of this integral 

G'±(C/,t) = (det(lTtC7fl))Ti (A.7) 

serves as a generating function for the characters of (anti)symmetrized products of arbitrary 
copies of the representation R. More explicitly, by expanding the determinant we may 
express the results directly in terms of the character for the representation R as 



oo 

n=0 
oo 



n=0 



Appendix B. Counting states in U{N) gauge theories 

In this appendix we derive the precise formula ( p.6|) for the counting of gauge-invariant 
states in a large theory with adjoint fields, and we discuss the single-particle partition 
functions for theories on a sphere with various field contents. 

B.l. Counting gauge-invariant states precisely 

In order to count the number of independent operators corresponding to traces of prod- 
ucts of fields in the large N limit, we wish to count the number of different arrangements 
of objects subject to an identification of arrangements related by a cyclic permutation. 
This can be done using Polya's theorem 

Consider a set of m types of objects, and associate a weight Xi with each of these 
objects. The weight associated with a collection of these objects is simply the product 
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of the weights associated with each of the individual objects. Polya solved the general 
problem of summing over weights for all sets of k of these objects, two sets being treated 
as identical if they are related to each other under the action of a specified subgroup of the 
permutation group. The subgroup relevant to us is simply the cyclic subgroup of order k; 
we will state Polya's result for this case. Define the polynomial 

Pk{yi,y2, ■ ■ ■ ,yk) - 2/2 ■■■Vk ' \^-^) 

TT 

where the summation in ( p.l| ) runs over all elements n of the cyclic subgroup, and n{7r)i 
is the number of cycles of length i in the permutation tt. The answer to the question 
addressed earlier in this paragraph is simply 

mm m 
1=1 1=1 i=l 

Applying this result to our problem, we find that the large partition function of 
single-trace states with k oscillators is precisely given by 

Zk = Pk{z{x), ■ ■ ■ , (B.3) 

where, as in §3, z{x) is the single-particle partition function. This implies that 

oo 

ZsT = ^Pk{z{x),z{x^),z{x^), ■ ■ ■ , z{x^)) = 

k (B-4) 

oo ^ fe ^ ' 



k 

k=l 1=1 

where n{k,l)q refers to the number of cycles of length q in the cyclic permutation by / 
shifts of k objects. 

It is easy to convince oneself that for specific values of k and Z, n{k,l)q is non-zero for 
only one value of q. At that value of q it is given by G{1, k), the greatest common divisor of 
/ and k. It then follows that the q for which l)q is non-zero is given hy q = k/G{l, k). 
Consequently, ( [B.3| ) may be rewritten as 

^ST = Y.\i2 ;2(x^/^(^''=))^(^'^). (B.5) 



k 

k=i 1=1 
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We now group together all terms with the same (fixed) q = k/G{l, k), so that k = G{1, k)q. 
Denoting j = G{l,k), we change the sum over / and /c to a sum over j and q, where 
each term appears once for every / < jq such that G{l,jq) = j. The number of such 
/'s is precisely (p{q), the number of positive integers which are not larger than q and are 
relatively prime to q. Thus, we obtain 

j = lg = l -"^ q=l 

as in (|3.6| ). 

B.2. Evaluating single-particle partition functions on spheres 

Next, we turn to a different topic which is the evaluation of the single-particle partition 
functions for d-dimensional field theories compactified on S'^~^ x IR, with unit radius for 
^d-i This may be carried out directly by noting that the Laplacian on the sphere (or 
the spatial parts of the other wave operators corresponding to particles with spin) may be 
written directly in terms of angular momentum generators, which may be diagonalized in 
the usual way. 

Alternatively, since free field theories are conformally invariant, and we are interested 
in conformally coupled fields (though it is easy to generalize our results also to other cases), 
we can use the conformal transformation that relates S'^~^ x IR to IR'^. This transformation 
takes states of the field theory on S'^~^ x IR to local operators on IR'^, with the energy 
of the state becoming the scaling dimension of the operator. Thus, an equivalent way to 
define the partition function in such a case is by z{x) = Yliiocai operators •'^^ where A is 
the scaling dimension of the operator. 

We begin by considering a free scalar field ip. The local operators in the theory are 99, 
diip, didjip, and so on, modulo the equation of motion. Ignoring the equation of motion 
for a moment, these operators are all generated by repeated application of the d different 
derivative operators 81,82 ■■ ■ 8d, each of which is of unit dimension and so has the partition 
function ji^^, on the free field if of dimension (d/2 — 1). Multiplying the various partition 
functions we find 

™d/2-l 

In order to obtain zs{x) we must subtract from this the partition function for operators 
that vanish by the equation of motion 8"^^ = 0. Such operators are generated by acting 



78 



with an arbitrary number of derivatives on 9^9?, so their partition function is x'^z'g(x). 
Thus, we find 

zs{x) = (1 - x^)z's{x) = (B.8) 

As a check, we note that in (i = 4 the operators that we get by acting with k derivatives are 
in the k^^ traceless symmetric representation of SO [A) which has ji = j2 = f ? and they 
have dimension A = /c + 1. These are simply the scalar spherical harmonics on . This 
implies that the number of operators of dimension A is ns(A) = (2ji + l)(2j2 + 1) = 
consistent with the Taylor expansion of ( [B.8[ ) for d = A. 



Next, we turn to the free vector field. The number of gauge-invariant operators is 
independent of the gauge, so we can fix an arbitrary gauge for the counting. We will use 
the gauge = on S'^~^ x IR, which becomes the gauge x^A^ = after the conformal 
transformation to IR*^ (recall that, according to the state-operator map, all operators are 
to be evaluated at a; = 0). Differentiating the gauge condition at the point a; = we find 
the relations 

= 0, d^A, + d^A^ = 0, ■ ■ ■ , d{,^d,, . . . A,^^^} = 0, ■ ■ ■ (B.9) 

where the brackets {} denote symmetrization. To start with we ignore both ( |B.9| ) and the 
equation of motion - this leads to a single-particle partition function z'y{x) = x'^~ i dz'g{x) 
(since the gauge field must have scaling dimension one in any space-time dimension). 
Operators of dimension A that are set to zero by (|B.9|) are given by symmetric ten- 
sors of rank A; based on the previous paragraph the corresponding partition function is 
x^~i z'g{x) — 1, where the last subtraction comes because there are no tensors of rank zero 
in (|B.9| ). With the condition ( [B.9| ), the Maxwell equation (at x = 0) simply reduces to 
d^Af^ = 0. The number of independent operators set to zero by the equation of motion 
is, therefore, counted by dx'^~ i z'g{x) . Finally, the number of operators set to zero by 
both the constraint ( [B.9| ) and the equation of motion is encoded in the partition function 
z'g{x). Putting it all together, using (|B.7|) , we find 

dx 1 dx^ x^ 

^y^"") = (i-xy ~ {i-xy " (i-x)'^ (T^^ ^ , , 

^ ^ ^ ' ^ ' ^ ' B.IO 
{l + x){l+x^ - dx) 

(l-x)^-i ■ 

As a check, we note that in four dimensions, the set of operators formed by acting with 
k derivatives on A^, obeying (|B.9| ) and the equation of motion, transform in the SO (4) 
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representation (ji, ^'2) = (^2^7 ■^^)©(^y^7 ■^2^)- These are the vector spherical harmonics 
on S^. It foUows that the number of operators at dimension A is ny(A) = 2(A^ — 1), 
consistent with (|B.1C1|) . 



FinaUy, we turn to free fermions. For concreteness we work in even dimensions with 
complex spinors of no chirality restrictions. Such a spinor has 2^+^ real components. Ig- 
noring the equation of motion, the partition function for spinors is z'p{x) = 2i~^^ ^/xz'g{x). 
The partition function that counts the operators which are set to zero by the Dirac equa- 
tion, is 2i~^^xi z'g{x). Subtracting the second from the first we find 



2 2' T. 2 



Of course, ( [13.11 ) should be divided by two for chiral spinors or real spinors, and by four 



for spinors that are both chiral and real. As a check on ( [B.llD , note that, in d = 4, the 



operators made from a complex chiral fermion field, at dimension /c + ^, transform in the 



SO{A) representation with (ji,j2) = (f 5 ■^1^); there are 2k{k + 1) such operators (the 
factor of 2 is because the spinors are complex), in agreement with (|B.11|) . 

Note that each of (gj), (|RTo|) , and ( pUD tends — s> 1 (the high temperature 
limit) to 

2J\fdof 
{1-xY 

where M'^"^ is the number of physical real degrees of freedom in the corresponding field. 
The formulas in this section, used for d = 4, imply the formula (|3.20| ) of §3. 



4x) ^ 7^— WIT' (B-12) 



Appendix C. Hagedorn transitions at finite string coupling 

There have long been speculations that at finite values of the string coupling Qs, one 
could have a phase transition at (or near) the Hagedorn temperature, with a different 
description of the theory at high temperatures. In this appendix we analyze, following 
p^ , the effect of turning on a small string coupling Qs on the partition function. Recall 



that, as reviewed in §2.1, the Euclidean partition function includes a winding mode W that 
becomes tachyonic above the Hagedorn temperature. On general grounds, the perturbative 
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effective action that describes the interaction between W and aU other modes (let us caU 
them (fn) takes the form 

s^ff = I d-x\ \dw\^ + - + V i-id^nY + + ■■■] + 



(C.l) 

^ / X ^ _ , ,n o^, , A \ 



If aU other fields (f>n are massive we can integrate them out, and remain with an effective 
action for W near the Hagedorn temperature. After rescaling W this takes the form 

= 1 / A- (law-p + - Dii^p + + ...)+.... (c,2) 

where 

7-2 

b = e (C.3) 

The dynamics of this theory depends crucially on the sign of b. For 6 > 0, Seff is the 
Landau Ginzburg free energy for a system that undergoes a second order phase transition 
at T = Th- In this case, for T < Th, Seff is minimized at = 0. The 'saddle point' 
contribution to the free energy vanishes and the leading 0{g^) contribution to ln(Z) is 
given by the free string theory partition function. However, at temperatures just above 
Th, Seff is minimized at 

\W\^= fC4) 



and the 0{gg ^) saddle point contribution to the free energy is 

V gl ■ 327r^a'^Tfjb' 



(C.5) 



On the other hand, for 6 < 0, Seff is the Landau Ginzburg free energy for a system 
that potentially undergoes a first order phase transition at a temperature lower than Th- 
The 0{gg) free energy below the phase transition temperature is once again computed 
by free string theory. In this case the high temperature behavior is dominated by large 
values of W, so it cannot be controlled in string perturbation theory; if there exists a high 
temperature saddle point, ^"[f is again of 0{g~'^), but its precise value depends on the 
details of the terms we denoted by "■ ■ ■" in (|G.2| ). 
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For string theory in flat space the massless dilaton is always one of the modes 99^, 
giving an infinite negative contribution to b in ( |C.3| ) , so one expects that if there is a phase 
transition in this case it would be of first order. Unfortunately, presumably in this case 
(and for any string theory with finite coupling in IR*^"^'^) no sensible high-temperature 
phase exists, both because there is a Jeans instability which cannot be ignored once the 
free energy is of order l/g"^ [0, and because we expect the high-energy spectrum to 
include black holes with a density of states growing faster than exponential (see, e.g., 
||871| ). However, while the formulas we wrote above strictly apply only to the flat space 
case, a similar Hagedorn behavior may be found in other spaces as well, including spaces 
like anti-de Sitter space where the spectrum effectively has a mass gap, and where these 
problems do not occur. We expect that a similar effective action would arise also in these 
casesi^l. In such cases, for example in type IIB string theory on AdS^ x S^, the arguments 
presented above may apply, with the value of b being either positive or negative, depending 
on the dynamics. It turns out that the specific case of string theory on anti-de Sitter space 
in global coordinates corresponds to the case of 6 < 0, since it is known to exhibit a 
first order phase transition which occurs (at small curvatures) well below the Hagedorn 



temperature [1T3|JT1[] . As expected from the discussion above, this phase transition is not 
visible in string perturbation theory around the AdS^ background but requires additional 
input. 

It is interesting to note the similarity between the analysis here and the analysis of §6 
of the Hagedorn transition in weakly coupled large gauge theories. 



Though, in theories that effectively have finite volume, we do not expect to find strict phase 
transitions at finite coupling. 
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